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Abstract 

We  propose  a  general  framework  for  obtaining  asymptotic  distributional  bounds  on  the 
stationary  backlog  WAl+A'2’c  in  a  buffer  fed  by  a  combined  fluid  process  A\  +  A2  and  drained  at 
a  constant  rate  c.  The  fluid  process  A\  is  an  (independent)  on-off  source  with  average  and  peak 
rates  p\  and  r i ,  respectively,  and  with  distribution  G  for  the  activity  periods.  The  fluid  process 
A-2  of  average  rate  p2  is  arbitrary  but  independent  of  A\.  These  bounds  are  used  to  identify 
subexponential  distributions  G  and  fairly  general  fluid  processes  A2  such  that  the  asymptotic 
equivalence  P  |lF  *'  ’  ’  >  x\  ~  P  \yvAl,c~P2  >  x\  (x  — >  oo)  holds  under  the  stability  condition 

pi+p-2  <  c  and  under  the  non-triviality  condition  c-p2  <  r\.  The  stationary  backlog  WAl,c~P2  in 
these  asymptotics  results  from  feeding  source  A\  into  a  buffer  drained  at  reduced  rate  c—p2.  This 
reduced  load  asymptotic  equivalence  extends  to  a  larger  class  of  distributions  G  a  result  obtained 
by  Jelenkovic  and  Lazar  [18]  in  the  case  when  G  belongs  to  the  class  of  regular  intermediate 
varying  distributions. 

‘This  work  was  performed  in  part  while  R.  Agrawal  was  visiting  INRIA  during  the  Spring  of  1997  and  while  A. 
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1  Introduction 


On-off  sources  provide  a  natural  and  versatile  tool  for  modeling  incoming  traffic  at  a  node  (router) 
or  gateway  of  a  network.  An  on-off  source  is  characterized  as  a  stochastic  process  which  alternates 
between  periods  of  silence  (off)  and  activity  (on).  During  on  periods  such  a  source  generates  data 
continuously  at  a  constant  rate,  and  becomes  silent  during  off  periods;  a  more  acurate  definition 
will  be  given  in  Section  2. 

As  traffic  flows  generated  by  multiple  on-off  sources  are  typically  multiplexed  onto  a  single  link,  it 
is  of  great  practical  importance  to  address  the  corresponding  buffering  issues  in  an  effort  to  make 
efficient  use  of  network  resources.  This  is  often  carried  out  in  the  context  of  the  following  simple 
model  for  a  multiplexer:  The  superposition  of  these  on-off  sources  is  offered  to  a  single  infinite 
capacity  buffer  which  is  drained  at  a  constant  rate.  If  W  denotes  the  resulting  stationary  backlog 
(assumed  to  exist),  then  it  is  expected  that  its  probability  distribution  crucially  depends  on  the 
statistics  of  the  on  periods.  For  instance,  it  is  well  known  that  for  a  single  exponential  on-off  source, 
i.e.,  a  source  with  on  (and  off)  periods  which  are  exponentially  distributed,  the  tail  distribution 
P  [W  >  x]  decays  exponentially  fast  as  x  tends  to  infinity  [1].  This  exponential  decay  property  is 
preserved  under  multiplexing  in  the  sense  that  when  several  independent  exponential  on-off  sources 
are  combined,  the  tail  distribution  P  [W  >  x]  still  decays  exponentially  fast  [12].  Both  situations 
discussed  so  far  are  instances  of  a  class  of  Markov  modulated  fluid  models  which  has  been  extensively 
studied  [1,  12,  23,  27]  since  the  seminal  work  of  Kosten  [20].  A  fairly  comprehensive  theory  has 
been  developed  for  such  sources,  and  algorithms  are  now  available  for  the  numerical  evaluation  of 
the  tail  distribution  P  [W  >  x]  for  all  values  of  x  [2,  27]. 

On  the  other  hand,  the  situation  is  quite  different  when  at  least  one  of  the  on-off  sources  has 
heavy-tailed  on  periods.  The  need  for  considering  such  models  with  heavy-tailed  components  can 
be  traced  back  to  recent  measurements  of  network  traffic  [21]  which  exhibit  long-range  dependence 
and  burstiness  over  an  extremely  wide  range  of  time  scales.  Along  these  lines,  for  a  single  on-off 
source  with  subexponential  distribution  G  for  the  on  periods,  Jelenkovic  and  Lazar  [18,  Thm  9]  have 
shown  that 

POO 

P  [W  >  x]  ~  K0  /  (1  -  G(u))  du  (x  ->  .30)  (1.1) 

for  appropriate  constants  A'o,«o  >  0  determined  by  the  source  statistics  (such  as  average  and  peak 
rates)  and  the  buffer  release  rate  (Remark  4.1);  typical  examples  of  subexponential  distributions 
include  the  Weibull,  log-normal  and  generalized  Pareto  distributions  [14].  In  fact,  the  asymptotic 
(1.1)  extended  a  result  obtained  earlier  by  Boxma  [6,  Thm  5.1]  for  generalized  Pareto  (or  regularly 
varying)  distributions  G  of  the  form 

1  —  G(x)  =  x  '/.(.(•!.  x  >  0 

with  0  >  0  and  some  regularly  varying  L  :  IR+  —  IR+  [5].  Recently,  the  case  1  <  0  <  2  has  been 
viewed  with  particular  interest  since  it  corresponds  to  the  input  process  being  long-range  dependent. 

Extensions  of  (1.1)  to  multiple  on-off  sources  have  been  considered  in  the  literature,  with  a  survey  of 
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related  recent  results  available  in  [7].  In  [6,  Thm  6.1]  Boxma  already  showed  that  if  a  single  on-off 
source  (source  1)  with  regularly  varying  on  periods  shares  the  buffer  with  an  exponential  on-off 
source  (source  2),  then  the  former  dominates  the  behavior  of  the  buffer.  This  result  was  extended 
by  Jelenkovic.  and  Lazar  [18,  Thm  10]  to  two  on-off  sources  with  more  general  statistics.  It  was 
shown  that  source  1  dominates  the  behavior  of  the  buffer  if  its  on  periods  have  an  intermediate 
regular  varying  distribution  [10]  as  long  as  the  tail  of  the  on  periods  of  source  1  is  heavier  (in  some 
precise  technical  sense)  than  that  of  source  2.  A  noteworthy  byproduct  of  the  results  in  [6,  18] 
is  that  the  source  2  contributes  to  the  asymptotic  behavior  of  P  [W  >  x]  only  through  its  average 
fluid  generation  rate  p2-  In  fact,  the  following  rephrasing  of  these  results  was  first  pointed  out  by 
Jelenkovic  and  Lazar  [18,  Thm  10]:  Let  A\  and  A 2  denote  the  on-off  sources  1  and  2  with  average 
rates  p\  and  p-j,  and  peak  rates  ?’i  and  i’2,  respectively.  The  combined  (or  multiplexed)  fluid  process 
A\  +  A-2  is  offered  to  a  buffer  which  is  drained  at  rate  c  (fluid  units/sec.)  under  the  stability  condition 
Pi  +  P-2  <  c.  Under  the  non- triviality  condition  c  —  p2  <  r  1,  the  corresponding  stationary  backlog 
WA 1+A2X  }ias  the  property 

p  [W-.4i+A2,c  >  _rj  _  p  [ir  '  >  ,r]  .  X  )  (1.2) 

where  the  stationary  backlog  IT  !  r  results  from  feeding  source  A\  into  a  buffer  drained  at  the 
reduced  rate  c  —  p2-  This  asymptotic  equivalence  (1.2)  reflects  the  following  intuitive  notion:  The 
tail  of  Pareto  distributed  activity  periods  of  source  1  is  considerably  heavier  than  the  exponential 
tails  governing  the  exponential  on-off  source  2.  As  a  result,  a  single  on  period  for  source  1  is  likely 
to  correspond  to  a  large  number  of  successive  on  and  off  periods  in  source  2.  Such  a  disparity  in 
time  scales  is  enough  for  the  Law  of  Large  Numbers  to  kick  in  for  source  2,  effectively  averaging  out 
random  fluctuations  about  the  mean  p2  and  replacing  them  by  the  average  behavior  of  the  source 
2.  It  is  now  a  small  step  to  believe  in  the  plausibility  of  (1.2). 

The  reduced  load  (asymptotic)  equivalence  (1.2)  suggests  a  natural  way  of  approximating  the  dis¬ 
tribution  of  WAl+A 2,c  with  that  of  WAl,c~p 2.  As  this  latter  quantity  is  associated  with  the  single 
source  T| ,  a  reduction  in  computational  efforts  may  result,  at  least  in  principle.  For  instance,  when 
applicable,  (1.1)  and  (1.2)  together  imply 

P  \WM +a2,c  >  ,x]  _  A' |  r  (l  _  G(u)}  du  (x  ->  00) 

^  Jx/OL  l 

for  appropriate  constants  A'i,ai  >  0  determined  by  the  statistics  of  A\  and  by  the  release  rate 
c  —  p2-  Given  its  asymptotic  basis,  this  approximation  will  become  increasingly  accurate  with  x 
large,  a  property  that  might  make  it  well  suited  in  various  contexts  for  evaluating  very  small  cell 
loss  probabilities  via  buffer  overflow  probabilities.  However,  the  accuracy  of  the  resulting  estimates 
remains  an  open  question,  with  some  indications  that  it  might  be  poor. 

Leaving  aside  these  computational  issues,  we  note  that  the  class  of  intermediate  regular  varying 
distributions  in  [18]  includes  regularly  varying  distributions,  but  does  not  contain  the  log-normal  and 
Weibull  distributions.  Hence,  the  validity  of  (1.2)  is  already  in  question  when  the  activity  period  of 
source  1  is  characterized  by  these  standard  subexponential  distributions.  However,  the  plausibility 
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argument  made  for  (1.2)  in  the  Pareto  case,  if  indeed  correct,  holds  out  the  possibility  that  the 
range  of  validity  for  (1.2)  extends  beyond  the  class  of  intermediate  regular  varying  distributions. 
In  particular,  in  the  same  way  that  Boxma’s  result  for  single  fluid  with  regularly  varying  activity 
periods  in  [6]  was  generalized  as  (1.1)  by  Jelenkovic  and  Lazar  [18],  it  is  natural  to  speculate  whether 
(1.2)  holds  more  generally  for  the  class  of  subexponential  distributed  activity  periods.  This  question 
forms  the  motivation  behind  the  developments  presented  here  as  we  revisit  the  model  in  [18]  when 
source  1  has  a  subexponential  activity  period. 


To  help  the  reader  navigate  the  many  technical  sections  of  the  paper,  we  summarize  below  the 
approach  taken  to  establishing  (1.2)  and  some  of  the  paper’s  main  contributions:  The  arguments 
articulate  around  lower  and  upper  bounds,  which  in  the  best  of  cases  take  the  form 


1  <  lirninf 


p  fir  1  ■ 


>  X 


and  lim  sup 


p  \wAl+A'2 


>  x 


Ull  __  r  A  Uiil  UULf  __ _ r  _  _  _  A  1 

P  [],]/AuC-p2  y  x ]  p  [p/U1,C-p2  >  >r] 

Our  point  of  departure  for  establishing  such  bounds  is  the  representation 

WAi+A2,c  t  sup(.l|(/)  +  A2(l)  -  ct) 

t>  0 


<  1. 


(1.3) 


(1.4) 


which  holds  under  the  usual  assumptions  (Section  2).  As  in  [18]  we  introduce  perturbations  h  : 
R+  — >  R  in  order  to  write  (1.4)  as 

Wai+a2x  =gt  sup (Ai{t)  -  (c-  p2)t  +  A2{t)  -  pit) 
t>  o 

=  sup  (Ai(t)  -  (c  -  p2)t  +  h(t)  +  A2(t)  -  pit  -  h(t)) .  (1.5) 

*>o 

While  only  linear  perurbations  sufficed  in  [18],  we  shall  need  general  perturbations  for  handling  the 
broader  class  of  sub  exponential  distributions.  This  decomposition  (1.5)  is  then  invoked  in  Section 
2  in  order  to  derive  generic  bounds.  These  bounding  arguments  hold  in  a  fairly  general  framework, 
and  are  given  in  terms  of  the  “perturbed"  backlog  associated  with  source  1  after  load  reduction, 
namely 

WM'C  :=  sup  (A1(t)  -(c-p2)t  +  h(t )) . 
t>o 

Asymptotic  bounds  on  the  tail  distribution  of  the  random  variable  (rv)  WAl,c~P2'h  are  established 
in  Section  4  when  the  source  A\  is  a  standard  independent  on-off  source  with  subexponential  ac¬ 
tivity  periods;  the  needed  facts  on  subexponential  distributions  are  collected  in  Section  3.  These 
asymptotic  bounds  on  WAl'c~p 2,h  can  now  be  used  in  conjunction  with  the  bounds  developed  for 
WA i+A#?  in  Section  2.  This  is  done  in  Sections  5,  6  and  7  when  the  source  A\  is  a  standard  inde¬ 
pendent  on-off  source  with  subexponential  activity  periods  and  source  A 2  is  a  fairly  general  fluid 
source  satisfying  at  minimum  a  Central  Limit  Theorem.  The  bounds  are  optimized  by  considering 
perturbations  of  the  form  eh  with  e  small  and  letting  s  go  to  zero  in  the  resulting  bounds. 

Still  we  need  to  identify  the  best  perturbation  “direction"  /i,  best  in  the  sense  of  making  the  lower 
(resp.  upper)  bound  largest  (resp.  smallest).  In  the  determination  of  a  suitable  h  for  either  the 
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R+  defined 


lower  or  upper  bound,  a  crucial  and  central  role  is  played  by  the  mapping  me  :  IR+ 
by 


ii)c(x) 


f~(l-G(t))dt 
1  -  G(x) 


x  >  0 


with  G  denoting  the  distribution  of  the  activity  period.  In  fact,  when  G  is  subexponential  but 
not  regularly  varying,  the  requirements  imposed  on  h  for  optimizing  the  bounds  coincide  with  well- 
known  properties  enjoyed  by  distributions  in  the  maximum  domain  of  attraction  of  the  Gumbel 
distribution  [14],  most  notably  properties  characterizing  the  so-called  auxiliary  function  associated 
with  such  distributions.  The  Weibull,  log-normal  and  Bentkander  distributions  belong  to  this  max¬ 
imum  domain  of  attraction,  thereby  allowing  us  to  take  li  (x )  ~  mg (x)  (x  — »  oo)  in  these  cases  and 
to  conclude  to  the  validity  of  (1.2)  under  certain  conditions;  these  examples  are  discussed  in  Section 
8.  The  importance  of  the  function  rnr;  is  reinforced  by  a  negative  result  of  Dumas  and  Simonian  [11] 
to  the  effect  that  the  asymptotic  equivalence  cannot  hold  if  liimE_00  mc(x)/ \fx  =  0.  In  particular, 
while  (1.2)  holds  for  all  log-normal  distributions,  it  holds  only  for  the  Weibull  distributions  that 
have  heavy  enough  tails  as  in  (8.17). 


Generalized  Pareto  rvs  belong  to  the  maximum  domain  of  attraction  of  the  Fechet  distribution,  and 
not  of  the  Gumbel  distribution  [14].  However,  the  approach  presented  here  applies  to  that  case  as 
well;  in  fact,  mc(x)  ~  x  ( x  — >  oo)  and  points  of  contact  with  Extreme  Value  Theory  also  emerge 
since  the  arguments  of  Jelenkovic  and  Lazar  [18]  were  based  on  using  linear  perturbations.  The 
somewhat  singular  nature  of  the  regularly  varying  case  is  made  increasingly  apparent  in  Section  9 
where  we  provide  various  extensions  for  it:  Following  [18]  we  show  that  the  validity  of  (1.2)  also 
holds  when  G  belongs  to  the  larger  class  of  intermediate  regular  varying  rvs  [9].  The  asymptotic 
“scale  invariance"  property  of  these  rvs  is  shown  to  imply  the  validity  of  the  upper  bound  in  (1.3) 
without  independence  of  the  sources. 

In  Section  10.1  we  consider  the  case  when  A\  is  a  superposition  of  independent  on-off  sources. 
Exact  asymptotics  for  this  case  are  notoriously  hard  to  obtain  in  general;  we  contend  ourselves  with 
a  lower  bound  that  extends  an  earlier  result  of  Choudhury  and  Whitt  [8].  Finally,  in  Section  10.2 
we  discuss  the  situation  often  encountered  in  practice  where  A-2  itself  is  an  aggregation  of  several 
sources,  typically  with  a  simpler  probabilistic  structure  (e.g.,  independent  on-off  sources);  we  show- 
how  the  necessary  technical  conditions  on  the  component  sources  transfer  to  the  aggregate  source. 
Section  11  closes  the  paper  with  a  list  of  open  problems  associated  with  the  reduced  load  equivalence 
(1.2). 


To  facilitate  the  reading  of  this  long  paper,  we  have  relegated  proofs  of  major  technical  results  to 
several  appendices.  A  word  on  the  notation  in  use:  Throughout  =>t  denotes  the  convergence  in  law 
with  t  going  to  infinity.  Equivalence  in  law  or  in  distribution  between  rvs  and  stochastic  processes 
is  denoted  by  =st ,  and  we  use  <st  for  the  strong  stochastic  ordering  between  rvs.  Also  for  mappings 

f,g  :  IR_|_  — *■  IR,  the  relation  f(x)  ~  g{x)  is  understood  as  lim.,. ^  =  1,  the  qualifier  (x  — *■  oo) 

being  omitted  for  the  sake  of  notational  simplicity.  For  any  scalar  x,  we  write  ,r+  =  max(,r,  0) 
and  with  any  mapping  tp  :  11+  — ►  IR,  we  associate  the  mapping  ip+  :  IR+  — ►  IR+  defined  through 
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y?+(a')  :=  (i^>(x))+  for  all  x  in  R. 


2  Generic  bounds 


We  adopt  the  following  framework  in  order  to  develop  several  generic  bounds:  A  fluid  process  (or 
source)  is  defined  as  any  R+-valued  stochastic  process  A  =  {A(t),  t  >  0}  with  non-decreasing 
right-continuous  sample  paths  such  that 


lim 

t — >oo 


m 

t 


=  p 


a.s.; 


as  usual  the  non-negative  constant  p  represents  the  average  rate  of  the  source. 

We  interpret  A(t)  as  the  amount  of  fluid  generated  by  the  source  in  the  interval  [0,  f).  If  the  fluid 
process  A  :=  {A(t),  t  >  0}  is  offered  to  an  infinite  capacity  buffer  from  which  it  is  drained  at  the 
constant  rate  of  c  (fluid  units /gee),  then  it  is  well  known  that  the  corresponding  backlog  at  time 
f  >  0  is  given  by 

WA,c{t)  =  sup  (A(t)  —  A(s)  —  c(t  —  .s)) 

o  <s<t 

provided  the  buffer  is  empty  at  time  f  =  0. 

The  following  facts  are  well  known  [22]:  If  the  arrival  process  A  has  stationary  increments,  then 
WA,c(t)  =>t  WA'C  where  the  stationary  backlog  rv  WA'C  has  the  representation 


WA'C  =  sup  (AR(t)  -  ct) 
t>  o 


(2.1) 


with  4/j  :=  {Aji(t),  t  >  0}  denoting  the  time-reversed  process  associated  with  A.  Moreover,  the  rv 
WA,C  is  a.s.  finite  under  the  stability  condition  p  <  c. 

The  fluid  models  encountered  in  practice  have  stationary  increments.  In  such  cases,  in  dealing 
with  (2.1),  with  a  slight  abuse  of  notation  we  denote  the  process  A^  by  A  instead,  or  equivalently, 
we  interpret  the  arrival  process  A  backwards  in  time.  We  note  that  in  many  important  instances 
(Section  4),  the  distributional  equivalence  {A^(t),  t  >  0}  =st  {-4(f),  t  >  0}  holds. 

Except  in  a  few  isolated  cases  (e.g.,  [1,  23,  27]),  characterizing  the  distribution  of  the  stationary 
backlog  rv  WA,C  is  a  difficult  task,  not  to  say  an  impossible  one.  As  a  result,  we  resort  to  studying 
the  tail  behavior  of  WA'C,  and  we  do  so  by  deriving  lower  and  upper  bounds  on  the  quantity  of 
interest  in  structured  situations.  In  particular,  we  have  in  mind  situations  where  several  fluid  sources 
are  multiplexed  onto  a  single  link.  As  the  bounds  are  obtained  by  a  perturbation  technique,  we 
find  it  useful  to  generalize  the  definition  (2.1)  by  associating  with  any  mapping  h  :  1R+  — *■  1R  the  rv 
WA'c'h  given  by 

WA'c'h  :=  sup  (-4(f)  -  ct  +  h(t))  . 

t>  0 
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In  due  course,  several  assumptions  will  be  imposed  on  such  a  perturbation  mapping  h. 

In  Sections  2.1  and  2.2,  we  consider  the  case  when  the  fluid  process  A  is  obtained  as  the  superposition 
of  two  fluid  processes  A\  and  A2  as  understood  earlier  in  this  section,  with  no  additional  assumptions. 
For  each  i  =  1,2,  let  pi  denote  the  average  rate  of  source  Ai  and  let  c  >  0  denote  the  release  rate 
of  the  fluid. 


2.1  Generic  lower  bounds 

We  begin  with  a  generic  lower  bound  on  the  tail  distribution  of  the  stationary  backlog. 


Proposition  2.1  Assume  the  arrival  process  A  to  be  the  superposition  A  =  Ai  +  A 2  of  two  inde¬ 
pendent  fluid  processes  such  that 

Aflt)  <  nt ,  t  >  0  (2.2) 

Then,  for  any  Borel  mapping  h  :  IR+  —  It,  it  holds  for  each  x  >  0  that 


W 


Ai+A2,c 


>  x\  >  inf  P 

\{t:(ri+p2-c)t+h(t)>x} 


[ A2(t )  >  p2t  +  h(t)\  j  P  |lF  '  r  >  x]  .  (2.3) 


Proof:  Fix  d  ^ >  1  and  pick  0  (5  1  —  d  Fix  x  ^  0.  There  exists  an  a.s.  finite  and  nonnegative 

rv  Ts  such  that 


mm 


in  (dx,  WAl,c  P2'hJ  <  (di(r*)  -  (c  -  p2)  t6  +  1i(ts))  . 


This  can  be  seen  by  considering  separately  the  cases  WAl  ,c  p2,h  _  ^  anc[  H  ’  ' :  •<  P2,h  < 

oo;  the  rv 

rb  may  depend  on  x. 


P  \w 


■  at  +  h(t) 

(t 

>  0)  with 

a 

:=  ri 

+ 

P2 

—  c.  The  rv  rb  being  independent  of  A 2,  it  holds  that 

A’c  >  x] 

PiAflrs) 

- 

ic¬ 

1 

to 

<3 

+ 

h{n) 

+ 

M 

>(rs)  -  p2T6  -  h{rs)  >  x] 

P[di(r6) 

- 

ic¬ 

1 

to 

+ 

Hts) 

> 

X, 

A2(t6)  -  p2r6  -  Ii(t6)  >  0] 

P[A-i(t6) 

- 

le- 

1 

to 

+ 

Hts) 

> 

X , 

a{n )  >  x,  A2(ts)  -  p2 rb  -  h(Tb)  >  0] 

P  [A2{t)  >  p2t  +  h{t)}  P  [Ai(t)  —  (c  —  p2)t  +  h(t)  >  x  |  Tf,  =  t]  P  [rb  £  dt] 


> 


> 


I {t :  a(t)>x} 


inf  P  [A2(t)  >  p2t  +  h(t)\  P  [di(r^)  -  (c  -  p2)rb  +  Ii(t6)  >  x] 

H:a(t)>x} 


inf  P  [Aflt)  >  p2t  +  h{t)\  P 

•>} 


inf  P  [A2{tf  >  p2t  +  h{t)}  P 


in  (clx,  W 


mm 


Ai,c—p2,h 


> 


1  -  6 


\yAlX-p2Jl  y 


X 


1-6 
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where  the  first  equality  made  use  of  the  constraint  (2.2)  on  source  A\,  and  the  last  inequality  follows 
from  the  definition  of  Tg.  Letting  8  j  0  in  this  last  inequality  yields  the  desired  conclusion.  I 

The  interest  in  this  lower  bound  resides  in  the  fact  that  the  two  sources  have  now  been  decoupled, 
with  source  A\  (resp.  A2)  entering  only  the  second  (resp.  first)  factor.  In  Section  4,  we  focus  on  the 
evaluation  of  the  first  factor  in  the  more  restricted  context  of  stationary  independent  on-off  sources. 


2.2  Generic  upper  bounds 

The  upper  bounds  derived  in  this  paper  all  flow  from  the  following  observation: 


Lemma  2.1  Assume  the  arrival  process  A  to  be  the  superposition  .4  =  A\+A2  of  two  fluid  processes. 
For  any  mapping  h  :  ]R+  —  R,  it  holds  that 

WA'C  <  ir 'i,r  <"'  h  +  wA'2'p2~h.  (2.4) 


Proof:  We  note  that 

WA'C  =  sup  (Ai(t)  -  (c  -  p2)t  +  h(t)  +  A2(t)  -  p2t  -  h(t)) 
t>  0 

<  sup  (Aft)  -  (c  -  p-2)t  +  h(t))  +  sup  (A-2(t)  -  p2t  -  h(t))  (2.5) 

t>  0  t>  0 


and  the  conclusion  (2.4)  follows.  I 

This  upper  bound  is  interesting  only  when  the  rvs  WAl'c~p 2,h  and  WA'2'P2'~h  are  a.s.  finite;  a 
necessary  condition  on  h  for  this  to  happen  is  given  by 


0  <  lirninf 

t — >00 


h(t) 


<  lirnsup 

t—>o O 


h(t) 


<  c-(pi+p2) 


with  the  condition  becoming  sufficient  if  both  outmost  inequalities  hold  as  strict  inequalities.  When 
some  of  these  outmost  inequalities  hold  only  as  equalities,  the  conclusion  depends  on  the  growth 
behavior  of  h  at  infinity,  such  situations  being  discussed  in  Section  6.3. 


3  Preliminaries 


For  easy  reference,  we  collect  below  some  definitions  and  technical  facts  that  are  used  throughout 
the  paper.  The  proofs  of  the  various  lemmas  are  provided  in  Appendix  A. 

We  first  recall  the  definitions  of  various  classes  of  probability  distributions  on  IR.+  which  are  of 
interest  here:  An  R+-valued  rv  A’  is  said  to  have 


a  long  tail ,  denoted  A'  G  C,  if 


lim 

x—>oo 


P  [A'  >  x  —  y] 
P  [A  >  ,r] 


1,  ye  1R; 


(3-1) 


a  subexponential  tail ,  denoted  A'  G  5,  if 


P  [A  +  A'  >  x] 
Inn  —  - - - 

x  ^°°  P  [A  >  x\ 


=  2 


where  A'7  is  an  independent  copy  of  A. 

Additional  material  on  these  classes  of  distributions  can  be  found  in  the  monograph  [14].  It  is  well 
known  that  the  class  5  is  a  subclass  of  C  [13].  Also,  it  is  worth  noting  that  (3.1)  holds  for  all  y>  0 
if  and  only  if  it  holds  for  all  V  <  o. 

We  follow  up  with  some  standard  (and  some  less  standard)  results  on  long-tailed  distributions. 


Lemma  3.1  Let  X ,  Y,  Z  and  T  denote  four  mutually  independent  1R -valued  rvs  such  that  X ,  Y 
and  Z  are  nonnegative. 

(1)  If  X  e  S  (resp.  C)  and  P  [Y  >  .r]  ~  cP  [A’  >  x]  for  some  positive  constant  c,  then  Y  G  S 
(resp.  C); 

(2)  If  X  G  C,  then  P  [A'  —  Y  +  d  >  x]  ~  P  [A'  >  x]  for  any  scalar  cl; 

(3)  If  Z  G  S ,  P  [A’  >  x]  ~  ci  P  [Z  >  .r]  and  P  [ T  >  x]  ~  c?  P  [Z  >  x]  for  constants  ci  >  0  and 
C'2  >  0,  then  P  [A'  +  T  >  x]  ~  (ci  +  C2)  P  [Z  >  x].  In  particular,  {X  +  T)+  G  5  if  ci  +  C2  >  0. 

Next,  with  any  !R+-valued  rv  A'  with  0  <  E  [A']  <  00,  we  associate  the  R+-valued  rv  A'*  whose 
distribution  is  the  integrated  tail  distribution  of  A',  namely 

P  [A'*  <  .)’]  =  f  P  [X  >  //]  du.  x  >  0.  (3.2) 

E  [A  J  J 0 

Some  useful  facts  concerning  A'*  are  contained  in  the  following  lemma. 


Lemma  3.2  Let  X  and  Y  be  independent  1R ^-valued  rvs.  If  X  G  C  with  0  <  E  [A]  <  00,  then  we 
have 


and 


P  [A'  >  x]  =  o  ( P  [A'*  >  x] ) 

P  [A  -  Y  >  u]  du  ~  E  [A]  P  [A*  >  a-]  . 


(3.3) 

(3.4) 
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The  next  two  lemmas  address  the  transfer  of  tail  properties  of  the  rv  A'  to  that  of  a  transformed  rv 
p(A')  for  some  mapping  p  :  1R+  — ►  1R+.  References  to  absolute  continuity  are  given  in  [19,  p.  336]. 

Lemma  3.3  Consider'  a  Borel  mapping  p  :  R+  —  1R.  such  that 

(i)  p  is  strictly  increasing  in  the  limit,  i.e.,  there  exists  xo  >  0  such  that  the  restriction  p  : 
[.ro,  oo)  — *■  R  is  strictly  increasing,  with  liin^oo  p(x)  =  oo  and 

<ho  :=  sup  <p(x)  <  oo;  (3.5) 

0<x<ixo 


(ii)  p  is  absolutely  continuous  on  (xo,oo)  and  the  limit 

lim  p'(x)  =:  $ 

x—>oo 


(3.6) 


exists  and  is  finite. 

Then,  for  any  1R  + -valued  rv  X  with  0  <  E  [A']  ,  E  [p+(A')]  <  oo,  we  have 

P[(p+(X)r>x]  E[A] 

P[p+(X*)>x]  ~  E  [<p+( A')] ' 


(3.7) 


Lemma  3.4  Let  p  :  R+  — >  ]R  be  a  mapping  strictly  increasing  and  convex  in  the  limit,  i.e.,  there 
exists  xo  >  0  such  that  the  restriction  p  :  [.ro,  oo)  — *■  R  is  strictly  increasing  and  convex.  Under  the 
finiteness  condition  (3.5),  it  holds  that  p(X)  £  C  (resp.  S)  if  X  £  C  (resp.  S). 


We  close  this  section  with  facts  that  will  help  us  identify  the  appropriate  perturbation  mappings 
needed  to  apply  the  generic  bounds  of  Sections  2.1  and  2.2.  For  any  ]R.+-valued  rv  A',  with 

p(X)  :=  sup  {x  >  0  :  P  [A'  <  x]  <  1}  , 

we  introduce  the  function  mx  ■  [0,//(A'))  — >■  ]R.+  given  by 

mx(x)  :=  E  [A]  >  y  >  x  G  I0’  h(X))-  (3-8) 

If  p(X)  =  oo,  then  p(X*)  =  oo  and  m\  is  defined  on  the  entirety  of  R+.  This  situation  is  not 

restrictive  for  our  purpose  as  we  have  fi(X)  =  oo  whenever  A’  £  C,  in  which  case  lim, _ mx(x)  = 

oo  by  Lemma  3.2. 


Lemma  3.5 


Consider  an  1R ^-valued  rv  X  with  p(X) 

lim  P  [A  >  x  +  yp(x)]  _ 
x-^oo  P  [A'  >  x\ 


-  oo  and  0  <  E  [A']  <  too.  Assume  that 
7  (y),  i/eR  (3.9) 
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(3.10) 

(3.11) 


Such  limits  are  invariant  under  asymptotic  equivalence,  i.e.,  if  tp  and  0  are  mappings  II  ,  — >  R , 
such  that  tp(: r)  ~  4’{x),  then  tp  satisfies  (3.9)  with  lim,.  .-x- /(.<•)  =  oo  if  and  only  i/>  does  with 
linilwco  tp(x)  =  oo ,  and  the  limiting  function  7  is  the  same. 

Limits  of  the  type  (3.9)  are  well  known  in  Extreme  Value  Theory  [14,  25]  where  they  occur  in  the 
characterization  of  maximum  domains  of  attraction;  we  refer  the  reader  to  [14,  Chap.  3]  [25,  Chap. 
1]  for  additional  information  on  this  topic.  Of  particular  interest  are  several  technical  facts  which 
are  summarized  below  for  easy  reference.  Recall  the  definition  of  the  Gumbel  distribution  A  as  the 
distribution  on  IR  given  by 

A(;r)  =  e~e  ,  x  G  R. 

The  needed  results  are  culled  from  [14,  Thm.  3.3.27,  p.  143],  [16,  Thm.  2.5.1],  [25,  Lem.  1.3,  p. 
41]  and  [25,  Cor.  1.7,  p.  46],  and  are  specialized  below  to  R+-valued  rvs  with  infinite  support. 


Lemma  3.6  The  R ^-valued  rv  X  with  //(A')  =  00  and  )  <  E  [V]  <  00  belongs  to  the  maximum 
domain  of  attraction  of  A,  denoted  X  G  MDA(  A),  if  and  only  if  there  exists  a  m.apping  tp  :  R  — ► 
(0,oo)  such  that  (3.9)  holds  with 


P  [A'  >  x  +  y<p{x )] 
P  [A  >  ,r] 


y  G  R. 


A  possible  choice  is  tp  =  mx  given  by  (3.8),  in  which  case  (3.11)  holds. 


The  log-normal,  Weibull  and  Benktander  distributions  belong  to  MDA(A)  [14,  pp.  149-150],  among 
others. 

In  view  of  Lemmas  3.5  and  3.6,  we  conclude  that  for  any  Revalued  rv  X  with  //(A')  =  00  and 
0  <  E  [A']  <  00,  membership  in  MDA{  A)  implies 


As  we  shall  see  in  Section  8,  the  conclusions  of  Lemma  3.6  hold  mutatis  mutandis  for  the  class  of 
generalized  Pareto  distributions;  interestingly  enough  these  are  exactly  the  rvs  which  are  in  the 
maximum  domain  of  attraction  of  Frechet  distributions  (0  >  0)  [14,  p.  121]. 


11 


4  On-off  sources 


4.1  Preliminaries 

An  independent  071-off  source  with  peak  rate  r  is  characterized  by  a  succession  of  cycles,  each 
such  cycle  comprising  an  off-period  followed  by  an  on-period.  During  the  on-periods  the  source  is 
active  and  produces  fluid  at  constant  rate  r  (unit  fluid/unit  time);  the  source  is  silent  during  the  off- 
periods.  The  on-period  durations  {Bn,  n  =  0, 1, . . .}  and  the  off-period  durations  {!„,  n  =  0,1,...} 
are  mutually  independent  sequences,  each  composed  of  i.i.d.  rvs  such  that 

0  <  E  [Bn\  ,  E  [In\  <  oo,  n  =  0,1,...  (4.1) 

It  is  convenient  to  introduce  the  sequence  of  epochs  { Tn ,  n  =  0,1,  .  ..}  marking  the  beginning  of 
successive  cycles,  namely  Tq  :=  0  and  Tn+ 1  :=  If.  +  Bj,  for  each  n  =  0,  l,.. ..  Thus,  at  time 

Tn  begins  the  (n  +  l)s*  cycle  with  of-period  duration  In  and  on-period  duration  Bn.  The  activity 
of  the  source  is  characterized  by  the  {0,  l}-valued  process  £  :=  {£(1),  t  >  0}  given  by 

OO 

m  :=  Y,  1  \Tn  +In<t<  Tn+ 1]  ,  1  >  0,  (4.2) 

n— 0 


with  the  source  active  (resp.  silent)  at  time  t  if  £(1)  =  1  (resp.  £(1)  =  0).  The  total  amount  of  fluid 
generated  in  [0,  t)  by  the  on-off  source  is  now  simply  given  by 

.4(1)  =  r  [  £(s)ds,  t  >  0.  (4.3) 

J  o 

The  integrability  condition  (4.1)  ensures  these  auxiliary  quantities  to  be  well  defined  and  finite. 
Both  processes  {-4(1),  t  >  0}  and  {£(!),  t  >  0}  have  right-continuous  sample  paths. 

Under  these  assumptions,  the  following  facts  are  well  known:  The  process  £  admits  a  (time)  sta¬ 
tionary  version  which  we  still  denote  by  {£(1),  t  £  IR}  (with  a  slight  abuse  of  notation).  Moreover, 
its  time-reversed  version  {£(— t),  t  £  11}  is  statistically  indistinguishable  from  {£(1),  t  £  ]R}  itself. 
Consequently,  the  distributional  equivalence  {Ajf(t),  t  >  0}  =st  {-4(1),  t  >  0}  does  hold  here  and 
(2.1)  takes  the  simpler  form 

WA’C  =st  sup  (.4(1)  —  ct)  (4.4) 

t>  o 

where  .4  =  {.4(1),  1  >  0}  computed  through  (4.3)  with  the  stationary  version  of  (4.2).  From  now 
on,  with  a  slight  abuse  of  terminology  we  refer  to  .4  so  defined  as  a  stationary  (independent)  on-off 
source. 

We  shall  find  it  handy  in  the  sequel  to  use  the  following  construction  of  this  stationary  (indepen¬ 
dent)  on-off  source:  We  postulate  rvs  {In,Bn,  n  =  0,1,...}  describing  the  alternating  sequence 
of  off-  and  on-period  durations  starting  with  an  off-period  of  duration  Iq;  if  Iq  =  0  the  source  is 
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construed  as  starting  in  an  on-period.  In  the  stationary  regime  considered  here,  standard  renewal- 
theoretic  considerations  require  that  (i)  (Io,  -Bo),  { /„ .  n  =  1, . . .}  and  {. Bn ,  n  =  1, . . .}  be  mutually 
independent  families  of  rvs;  (ii)  the  rvs  {Ini  n  1. . . .)  (resp.  { Bn .  n  =  1, . . .})  be  i.i.d.  rvs  with 
I\  (resp.  B\)  distributed  as  the  generic  off-period  (resp.  on-period),  and  (iii)  the  relations 

l(Io,B0)  |  Io  >  0]  =st  (/|  •  B| )  and  [(Io,-Bo)  |  Io  =  0]  =st  (0 ,B*)  (4.5) 

hold  with  If  independent  of  B\.  Under  such  assumptions,  we  check  that 

[Io  +  -Bo  |  Iq  >  0]  =st  If  +  B\  and  [Iq  +  -Bq  |  Iq  =  0]  =<+  -Bf  (4.6) 


and 


p  [Jo  =  o]  =  „  ,  „E,  |Clf  , ,  1  --P- 


The  average  rate  p  of  the  source  is 


E  [Bx]  +  E  [h 


y  A(t) 

p  :=  lim  - =  pr 

t—*o O  t 


and  an  independent  on-off  source  is  a  particular  instance  of  a  fluid  process  as  defined  in  Section  2. 


The  construction  developed  here  for  a  stationary  on-off  source  differs  from  the  usual  one  (e.g.,  [17]) 
in  two  respects:  The  first  cycle  always  starts  at  time  t  =  0  with  an  off-period  (albeit  of  possible 
duration  Io  =  0)  so  that  every  cycle  contains  an  activity  period.  Moreover,  the  model  is  prescribed 
through  the  requirement  (4.5),  instead  of  the  more  usual  requirement  (4.6)  (implied  by  it).  These 
features  will  simplify  the  presentation  and  discussion  of  several  results  by  permitting  direct  sample 
path  arguments,  notably  in  Proposition  4.1. 


4.2  Bounds  for  on-off  sources 


As  is  already  apparent  from  the  generic  bounds  discussed  in  Sections  2.1  and  2.2,  it  will  be  useful 
to  investigate  the  tail  behavior  of  WA,1,h  under  a  wide  class  of  perturbations  h  :  11+  — »•  11.  This 
will  be  done  by  first  establishing  bounds  on  WA'1'h,  and  then  by  identifying  the  tail  asymptotics  of 
the  bounds. 

With  any  mapping  h  :  11+  —  11,  we  associate  the  auxiliary  mappings  a,  b  :  11+  —  1R.  given  by 

</(/):  (r  •r: )/  •  //(/)  and  b(t)  :=  — yf  +  h(t) t  t>  0.  (4.7) 

The  R- valued  rvs  k  =  0, 1, . . .}  can  now  be  defined  by 


Xk  :=  a(Bk)  +  6(4),  1  =  0,1,... 


and  set 


y-4+,/j 


max 


h(  0), 


sup 

n  — 0.1.. 


k= 0 


A, 


(4.8) 


(4.9) 


13 


Under  the  enforced  assumptions,  the  rv  A'o  is  independent  of  the  i.i.d.  rvs  {A'^  k  =  1,  2, . . .}. 

A  mapping  h  :  It+  —  It  is  said  to  satisfy  condition  (Hr),  r  =  1,  2,  if 

(HI)  h  is  absolutely  continuous  on  It+  with  derivative  h'  satisfying  7  —  r  <  h'(t)  a.e.  on  IF!  y 
(H2)  h  is  absolutely  continuous  on  ]R.+  with  derivative  h'  satisfying  h'(t)  <  7  a.e.  on  It+. 

Proposition  4.1  Consider  a  Borel  mapping  h  :  It+  —  It. 

1.  If  h  is  superadditive ,  then  VA,1,h  <  WAn'h; 

2.  If  both  (HI)  and  (H2)  hold,  and  h  is  subadditive,  then  W A'1'h  <  VA'1,h . 

A  proof  of  Proposition  4.1  is  available  in  Appendix  B. 

4.3  Tail  asymptotics  of  the  rv  VAr,J' 

With  the  help  of  Proposition  4.1,  we  can  bound  the  tail  of  the  rv  WA'1'h  with  that  of  VAr,'h.  Hence, 
of  particular  interest  are  the  following  asymptotics  of  the  tail  distribution  of  the  rv  V A,'i  h . 

Proposition  4.2  Consider  a  mapping  h  :  It+  —  ]R.  which  satisfies  (H2)  together  with  the  conditions 

(H3)  The  limit  lim^^oo  h'(x)  =:  h'(oo)  exists  and  is  finite; 

(Hf)  The  mapping  a  given  by  (4-7)  is  strictly  increasing  in  the  limit  with  linx^oo  a(x)  =  00. 
Further ,  assume  the  following  conditions: 

(H5)  The  rv  A'i  is  integrable  with  E  [A’i]  <  0/ 

(H6)  a+(Hi)  G  C; 

(H7)  a+(B()  G  S. 

Then,  it  holds  that  VA'1'h  G  S  with 

P  [vAl'h  >x\~  (p  +  K(h))P  [a4)^)  >  x]  (4.10) 

where 

K(h)  :=  (( r  -  7)  +  A'(oo))-UAL.  (4.11) 
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A  proof  of  this  result  is  given  in  Appendix  C.  Two  useful  consequences  emerge  by  combining  the 
bounds  of  Section  4.2  with  the  asymptotics  of  Proposition  4.2.  Indeed,  Propositions  4.1(1)  and  4.2 
together  imply  the  following  lower  bound  asymptotics  on  the  tail  of  WA'h'1 . 


Corollary  4.1  For  any  superadditive  mapping  h  :  1R.+  —  ]R.  which  satisfies  (H2)-(H7),  it  holds  that 


lim  inf 

x—>oo 


p  wA'~>'h  >  x 
P[a+{B$)  >  x] 


>  p  +  K(h). 


Upper  bound  asymptotics  on  the  tail  of  W  l  /'  '  can  also  be  obtained  once  we  combine  Propositions 
4.1(2)  and  4.2.  We  omit  a  formal  statement  of  this  result  as  we  do  not  use  it  in  that  form  in  the 
sequel. 


4.4  Comments 

We  close  this  section  with  remarks  that  will  be  useful  in  the  sequel: 

Remark  4.1  We  note  that  h  =  0  satisfies  the  assumptions  of  both  parts  of  Proposition  4.1,  whence 
WAri  =  l  ’  !  ".  Moreover,  the  assumptions  of  Proposition  4.2  are  automatically  satisfied  provided 
B\  G  C  and  B*  €  S,  in  which  case  specializing  (4.10)  yields  the  exact  asymptotics 

P  111'  1  '  >  x]  ~  (1  ~p^pp  [(?’  -  7)R*  >  x]  .  (4.12) 

L  J  7  -  P 

This  result  was  first  obtained  by  Jelenkovic  and  Lazar  [18,  Theorem  9]  by  resorting  to  the  Palm 
theory  of  stationary  processes.  This  is  to  be  contrasted  with  the  direct  approach  taken  here. 


Remark  4.2  In  the  assumptions  of  Proposition  4.2,  if  we  add  in  (H4)  the  requirement  that  the 
mapping  a  is  convex  in  the  limit,  then  (H6)  and  (H7)  are  implied  by  the  conditions  (H6bis)  and 
(H7bis),  respectively,  with 

(H6bis)  B\  G  £; 

(H7bis)  B\  G  5. 

This  is  a  simple  consequence  of  Lemma  3.4  applied  to  p  =  a+ . 


Remark  4.3  A  convex  (resp.  concave)  mapping  h  :  R+  — >  1R.  with  h( 0)  =  0  is  absolutely  contin¬ 
uous,  and  necessarily  superadditive  (resp.  subadditive).  Consequently,  a  convex  (resp.  concave) 
perturbation  function  h  is  a  natural  choice  when  considering  the  lower  (resp.  upper)  bound.  In 
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fact,  many  of  the  needed  conditions  appearing  in  Propositions  4.1  and  4.2  are  then  easily  verified 
with  the  help  of  Jensen’s  inequality  together  with  the  observation  that  the  limits 


lim 

t — 


h(t) 

t 


lim  h'(t)  =  h'(oo) 


take  place  monotonically  (thus  exist)  and  are  finite.  Such  a  discussion  is  provided  in  Section  8. 


5  Towards  bounds  for  two  independent  sources 


We  now  deal  with  the  situation  where  the  arrival  process  .4  is  the  superposition  of  two  independent 
fluid  processes  .4i  and  .  1-2,  i.e.,  A  =  Ai  +  .42,  where  .4i  is  a  stationary  “independent  on-off"  as 
understood  in  Section  4.1  and  .42  is  arbitrary.  Source  .4i  has  peak  rate  ri  >  c  —  p-2  and  its  generic 
activity  period  B\  has  the  property  that  B\  G  C  and  B\  G  S.  Its  generic  inactivity  period  I\  has 
an  arbitrary  distribution  with  the  only  requirement  that  0  <  E  [Ii]  <  00. 

Before  presenting  the  lower  and  upper  bounds  in  Sections  6  and  7,  we  pause  to  introduce  some 
notation  that  simplifies  the  presentation  of  the  results:  In  the  context  of  the  two  sources  described 
above,  we  write 

a  :=  n  +  p-2  -  c  (5.1) 

and  with  any  mapping  h  :  IR+  — *■  IR,  we  recast  the  earlier  definitions  (4.7),  (4.8)  and  (4.11).  The 
mapping  a  :  IR+  —  IR.  is  given  by 


a(x)  :=  ax  +  h(x),  x  >  0, 

(5.2) 

and  we  set 

K(h)  :=(a  +  h' (°o))-^piL 

(5.3) 

with 

A'i  :=  aB!  -  (c  -  p2)h  +  h(Bi)  +  h(h). 

(5.4) 

We  also  write 

L{h):-\iP-p)PP12){1>  +  ,m) 

(5.5) 

with  p  :=  E  [B\\ 

(E  [i?i]  +E  [Zi] ) — 1  denoting  the  stationary  probability  that  source  1  is  active. 

Note 

that  Z(0)  =  1  if  h  =  0. 


Let 


R(x ;  h) 


P  [a  +  TO  >  X\ 
P  <  1  /if  >  x] 


x  >  0. 


Under  (H4),  we  see  that 

R-(h)  :=  lim  inf  R(x\  h) 

X^OO 


p kWWl 

—00  P  [aB$>a+(x)] 
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(5.6) 


=  lim  inf 


P  [B\  >  x] 


and  similarly, 


R+(h)  :=  lim  sup  R{x\  h)  =  lim  sup 


x^oc,  p  [aB\  >  a{x)} 
P  [Bt  >  x] 


;r>'  x  P  (1  /i,  >  d{x)\ 


(5.7) 


We  draw  the  reader’s  attention  to  the  similarity  of  the  limits  (5.6)  and  (5.7)  with  those  guaranteed 
by  Lemmas  3.5  and  3.6  (with  particular  reference  to  (3.12)).  This  connection  will  be  exploited  to 
identify  the  appropriate  perturbation  functions  h  in  the  bounds  developed  thus  far. 

While  the  quantities  L(h),  R-(h)  and  R+(h)  will  help  quantify  the  impact  of  the  first  source  Ai,  the 
generic  lower  bound  (2.3)  suggests  that  the  contribution  of  the  source  A2  will  be  expressed  through 
the  quantity 

(5.8) 


A (h)  :=  lim  inf  P  [^(f)  >  p2i  +  h(t )]  . 

t—>o O 


Its  alternate  expression 


A (h)  =  lim  inf  P 


t—>o O 


Mi)  h(t) 

~  P  -2  > 


t 


t 


indicates  already  the  possibility  that  its  value  will  be  determined  by  refinements  to  the  assumed 
“Law  of  Large  Numbers" 


n.  Mi) 

Inn  -  =  p2  a.s. 

t— *•  00  t 


(5.9) 


defining  the  average  rate  P2  for  source  A 2.  Such  refinements  include  the  Central  Limit  Theorem 
(CLT)  in  the  form 

'Mi) 


Mi)  ^ 


t 


p  2 


all 


(5.10) 


with  a  >  0  and  U  denoting  a  Gaussian  rv  with  zero  mean  and  unit  variance.  This  condition  is 
not  prohibitive  for  it  holds  in  great  generality  for  a  variety  of  on-off  sources  (as  implied  by  similar 
results  on  renewal  processes  [15])  and  for  superpositions  thereof.  Moreover,  the  relation 


P[A2(t)-p2t>h(t)]=P 


Z2(t )  > 


Hi) 

Vi 


t  >  0 


(5.11) 


points  to  the  need  to  impose  constraints  on  the  behavior  of  for  large  t  in  order  to  get  non-trivial 
limits  in  (5.8). 


6  Upper  bounds  for  two  independent  sources 

6.1  A  basic  upper  bound  for  two  sources 

We  begin  with  an  intermediate  result. 
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]R  which  satisfies  (Hl)-(H 7)  [with 


Proposition  6.1  Consider  a  subadditive  mapping  h  :  It+ 
7  <-e-  P2  and  p  ^  pi].  If 


P  hr  <’  >  J  =0  (P  [a+(Bf)  >  x])  , 


(6.1) 


then  it  holds  that 

p  \wA  1+A2’c  >  x] 

lj77p  >  ]  -  L(h>R+ih> 

with  L(h )  and  R+(h)  given  by  (5.5)  and  (5.7),  respectively. 


(6.2) 


In  (6.1),  the  finiteness  of  the  rv  WA2,P2’  h  is  implicitly  assum.ed  for  it  is  not  necessarily  guaranteed 
under  the  conditions  (H1)-(H7)  imposed  on  h.  This  finiteness  issue  and  the  tail  behavior  of  the  rv 
WA^~h  are  quite  delicate  when 


h(t) 

Inn  — —  =  0 

t—»  OO  / 

as  occurs  in  many  interesting  instances.  This  point  is  explored  in  some  detail  in  Section  6.3. 


(6.3) 


Proof:  Combining  Lemma  2.1  and  the  upper  bound  of  Proposition  4.1(2)  we  find 

WA1+A2,c  <  yAuc-p2.h  _|_  ■^/■A2,p2-h  (g .4) 

where  the  rvs  VAl'c~P2jh  and  IT  h  are  taken  to  be  independent.  From  Proposition  4.2  [with 
•/  <  c  p-2  and  p  <—  pi],  it  holds  that  VA'J'h  G  S  with 

p  |r  '  /■"•''  >  ;r]  _  (p  +  K(h)) P  [a+(.B*)  >  x] 

where  K(h)  is  given  by  (5.3). 

Parts  (1)  and  (3)  of  Lemma  3.1  readily  ensure  under  (6.1)  that 

P  |p  *  •'  +  wM^-h  >  x]  ~  (p  +  K(h)) P  [a+(B[]  >  x]  , 


and  it  immediately  follows  from  (6.4)  that 


lim  sup  ■ 


WA  1+Aa 


>  x 


<P  +  K(h). 


;r - P  [(l^(H^)  >  x] 

Remark  4.1  applied  to  source  Ai  [with  7  <—  g  —  p-2  and  p  <—  pi]  yields 

(1  ~p)pi 


(6.5) 


P  \wAt 


C-p  2 


>  X 


-P  [aBf  >  x] 


c-{pi  ■  p->) 

and  the  desired  conclusion  (6.2)  follows  upon  combining  (6.5)  and  (6.6). 


(6.6) 
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6.2  An  improved  upper  bound  for  two  sources 


Proposition  6.1  begs  the  question  of  how  to  choose  the  perturbation  mapping  h.  In  view  of  the  form 
of  the  upper  bound  (6.2),  an  obvious  criterion  for  selecting  h  is  that  L(h)R+(h)  be  made  as  small  as 
possible  in  order  to  yield  the  best  upper  bound.  To  gain  some  insights  on  how  this  could  be  achieved, 
we  note  that  if  h  assumed  non-negative  values,  i.e.,  h  :  It+  — *■  IR.  .  then  amongst  its  scaled  versions 
{eh,  e  >  0},  the  smallest  value  of  K(eh)  is  achieved  for  e  =  0  and  that  L( 0)  =  limfj q  L(eh)  =  1. 
This  remark  leads  to  the  following  improvement  to  Proposition  6.1;  its  proof  is  straightforward  and 
therefore  omitted  for  the  sake  of  brevity. 


Proposition  6.2  Consider  a  subadditive  mapping  h  :  1R.+  —  It  such  that  whenever  0  <  e  <  e*  for 
some  e*  >  0,  the  scaled  mapping  eh  satisfies  (Hl)-(H 7)  [with  7  <—  c  —  p2  and  p  <—  p\[.  If 

P  [WA2,P2-Sh  >  =  0  (P  [af(Bf)  >  4sJ)  ,  0  <  e  <  e*  (6.7) 

with  the  mapping  as  :  IR  1  — ►  It.  given  by  a.  Ur)  :=  at  +  eh(x)  (x  >  0),  then  it  holds  that 


lim  sup 

X^OO 


p  \wAi+A2 


>  x 


P  [WAl’c~P2  >  x] 


<  lim R+(sh). 

s  JO 


(6.8) 


In  order  to  improve  the  upper  bound  (6.8)  we  need  only  select  the  perturbation  direction  h 
that  makes  liimjo  R+(fih)  as  small  as  possible.  For  h  >  0,  we  see  that  R+{eh)  >  1,  whence 
liimjo  R+(eh)  >  1,  and  it  is  therefore  tempting  to  seek  h  :  1R+  — *■  1R+  such  that  limfj0  R+(eh)  =  1, 
in  which  case  (6.8)  becomes 


lim  sup 

X^OO 


p  ir  '  •  ,2  r  >  x 

P  [WAl'c~P2  >  x] 


<  1. 


(6.9) 


But,  by  virtue  of  Lemmas  3.5  and  3.6  (and  (3.12)),  whenever  B 1  G  MDA( A),  it  holds  that 


R+(emBl) 


lim 

X — >00 


P  [Bf  >  x] 


P  [Bf  >  x  +  «  1emB1(x)\ 


=  e 


£  >  0 


and  taking  h  =  /»/;. .  we  indeed  get  lim s^0R+(eh)  =  1,  whence  (6.9),  provided  the  appropriate 
assumptions  are  satisfied.  Thus,  in  most  cases  of  interest,  niB1  (or  an  asymptotic  equivalent)  is 
expected  to  be  the  perturbation  function  of  choice  for  getting  the  best  possible  upper  bound  (6.9) 
in  Proposition  6.2. 


6.3  On  the  condition  (6.1) 

In  order  to  assess  the  range  of  applicability  of  Propositions  6.1  and  6.2,  we  need  to  focus  on  condition 
(6.1)  which  quantifies  the  situation  where  the  second  source  has  a  “lighter  tail"  than  source  A\.  To 
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that  end,  consider  a  mapping  h  :  ]R.+  —  ]R.  such  that  the  limit 

lim  =:  L  (6.10) 

t—*o O  t 

exists  and  is  finite.  The  corresponding  rv  WA'2'p2'~h  will  be  a.s.  finite  (resp.  infinite)  if  L  >  0 
(resp.  L  <  0,  a  case  of  no  interest  here).  However,  if  instead,  L  =  0  (or  equivalently,  condition 
(6.3)  holds),  then  there  is  no  a  priori  guarantee  that  the  rv  WA2,P 2~h  will  be  a.s.  finite  as  already 
indicated  by  the  following  result. 

Lemma  6.1  Assume  the  CLT  refinement  (5.10)  to  hold  for  source  A 2.  Then ,  for  any  mapping 
h  :  R_|_  —  1R.  such  that  the  limit 

lim  inf  —-1  =:  jj  (6.11) 

t-+0 o  \ft 

is  finite,  we  have  WA2'P2~h  =  00  a.s. 


If  h  does  satisfy  (6.10)  and  (6.11),  then  L  =  0  necessarily. 
Proof:  Fix  t  z’  0  and  x  in  H.  Wg  note  that 


Vt  (z2(t)  -  ^0  <  wA2'P2'-’\ 


whence 


<  P  [wA^~h  >  X\ft  . 


Letting  t  go  to  infinity  along  a  sequence  for  which  the  the  liminf  in  (6.11)  is  attained  in  this  last 
inequality  and  invoking  (5.10),  we  find 


P  [all  -  H  >  x]  <  lim  P  h  >  xy/t ]  =  P  \w 

t — >00  J 


A2'p2~h  =  00 


(6.12) 


Finally,  we  get  the  desired  conclusion  upon  letting  x  go  to  —too  in  (6.12).  I 

Consequently,  for  the  rv  WA2,P 2~h  to  be  a.s.  finite  under  (5.10),  it  is  necessary  that  H  =  00.  We 
now  turn  to  finding  sufficient  conditions  on  h  under  which  the  rv  WA2'P2~h  is  a.s.  finite.  In  the 
process,  we  identify  its  tail  behavior,  thereby  providing  the  means  to  check  (6.1). 

The  discussion  will  be  carried  out  in  the  following  regenerative  framework  which  contains  most  fluid 
models  discussed  in  the  literature,  including  the  independent  on-off  sources  of  Section  4:  Source 
A-2  is  a  fluid  process  {.42(f),  t  >  0}  process  characterized  by  a  succession  of  cycles,  where  for  each 
n  =  1,  2, . . .,  the  (n  +  l)st  cycle  has  duration  Cn  (with  Cn  >  0  a.s.)  and  is  associated  with  source  .42 
producing  fluid  in  amount  Yn  (with  Yn  >  0  a.s.).  Alternatively,  we  have  Yn  :=  A-2(Tn)  —  A-i^Tn-i) 
where  Tn  denotes  the  beginning  of  the  nth  cycle  (with  the  convention  Tq  =  0),  i.e.,  Tn  =  C\  +  . .  -+Cn. 
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No  additional  details  on  the  operation  of  source  A2  will  be  needed.  We  assume  that  the  rvs 
{(Yn,Cn),  n  =  1,2,...}  are  integrable  and  mutually  independent  with  {(!’„, Cn),  n  =  2,3,...} 
being  identically  distributed  rvs.  By  the  Renewal  Reward  Theorem,  the  mean  rate  p-j  given  by  (5.9) 
can  be  evaluated  as  P2  =  E  [I2]  / EfCy.  A  fluid  process/source  satisfying  the  above  requirements 
will  be  called  a  regenerative  fluid  process/source. 

The  other  key  probabilistic  assumption  is  the  existence  of  finite  exponential  moments:  There  exists 
a  constant  do  >  0  such  that 

Mn{6)  :=  E  [e0CVi]  <  .30  and  Nn{6)  :=  E  [e01"]  <  .30,  0  <  6  <  6»0,  n  =  1,  2.  (6.13) 


A  regenerative  fluid  process  satisfying  (6.13)  also  admits  the  CLT  refinement  (5.10). 


Proposition  6.3  Let  {.42(f),  t  >  0}  be  a  regenerative  fluid  process  satisfying  (6.13).  Consider  a 
nondecreasing  mapping  h  :  R+  — *■  ]R+  which  satisfies  both  (6.10)  with  L  =  0  and  (6.11)  with 
H  =  00.  Define  the  mapping  g  :  It+  —  R+  by 

g(t)  =  \ft  inf  — t  >  0. 

s>t  s 

(1)  Then,  there  exist  finite  constants  71,72  >  0  and  73  >  0  such  that 

P  I 'wA2,P2-h  >  J 

lim  sup  —^5 - - - -  <73-  (6-14) 

•E^°°  /  e“729  (t)/t  dt 

J-yix 


(2)  Furthermore,  if  there  exists  no  in  (0,1/2)  such  that 

r  •  t  h (*) 

Inn  inf  —7— —  =  00, 

t—  00  fi/2+i'o 


then  there  exist  finite  constants  6\  >  0  and  62  >  0  such  that 

p  [' [yA2,P2,-h  y  _rj 


lim  sup  ■ 


f.1  —  12  p  —  b\Xv 


<62,  v  G  (0,  2vo)- 


(6.15) 


A  proof  of  this  result  is  available  in  Appendix  D.  Proposition  6.3  provides  a  natural  vehicle  for 
checking  condition  (6.1)  as  is  done  in  Section  8  on  a  variety  of  examples.  Note  that  the  denominator 
in  (6.15)  is  vanishingly  small  with  x  large,  so  that  this  result  does  yield  a  non-trivial  bound  on  the 
tail  of  WA2-p2'~h. 
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The  situation  where  L  >  0  in  (6.10)  can  be  handled  in  a  variety  of  ways:  Indeed,  under  such  an 
asymptotic  linearity  assumption  on  h,  there  exist  constants  n/,  and  o'/,  >  0  such  that 

/'(/)  >  ah  +  / 3h.t ,  t  >  0 


and  the  comparison 


WA2'p2'~h  <st  WA2’P2+Ih  -  ah 


follows.  Whenever  the  source  A2  belongs  to  the  class  of  Markov  modulated  fluid  sources  [12]  or  more 
generally,  is  an  exponential  source,  as  understood  in  [18],  then  the  tail  behavior  of  WA2'p2~h  is  at 
most  exponential  since  that  of  WA2,p2+^h  has  exponential  decay.  On  the  other  hand,  if  A2  is  an  on-off 
source  with  subexponential  activity  periods,  then  the  asymptotics  (4.12)  [with  7  <—  P2  + 1%  ,  r  <—  /921- 
can  be  invoked.  In  either  case,  these  remarks  serve  as  the  basis  for  checking  (6.1). 


7  Lower  bounds  for  two  independent  sources 

7.1  A  basic  lower  bound  for  two  sources 

Proposition  7.1  Consider  a  superadditive  mapping  h  :  IR+  —  It  which  satisfies  (H2)-(H 7)  [with 
7  <—  c  —  p2  and  p  <—  pi].  Then,  it  holds  that 


L(h)A(h)R-(h)  <  lirninf 

x— >oo 


P  WA!+A2,c  y 

P  [WAl'c~P2  >  x] 


with  I.(li) ;  A (h)  and  R-(h)  given  by  (5.5),  (5.8)  and  (5.6),  respectively. 


(7.1) 


Proof.  Fix  x  >  0.  We  conclude  from  Proposition  2.1  applied  to  h  that 

p  [ir  1  •'  >  :r]  P  [ir  '  •  >  s] 


inf  P  [A2(t)  >  Pit  +  h(t)]  R{x:  h)- 


{t:a(t)>x} 


P  [a+(B()  >  x] 


P  [aB\  >  x] 


(7.2) 


As  we  have  in  mind  to  let  x  go  to  infinity  in  (7.2),  we  note  the  following:  By  Corollary  4.1  applied 
to  h  with  the  source  A\  [and  7  =  c  —  p%\,  it  holds  that 


p  +  I\(h )  <  lirninf 

X^-OO 


P  IP  1  r  P1,h  y  x 

P[a+(B[)  >  x] 


(7.3) 


with  K(h)  given  by  (5.3).  Moreover,  under  (H4),  the  mapping  a  is  eventually  strictly  increasing 
(thus  invertible)  with  Inn,.  .  N  a(x)  =  00 ,  and  we  readily  check  that 


lirninf  inf  P  [.42(f)  >  p2t  +  h(t)]  =  liminfP  [.42(f)  >  p2t  +  h(t)]  =  A (h).  (7.4) 

S-’^-oo  {t:a(t)>x}  00 
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Thus,  letting  x  go  to  infinity  in  (7.2)  and  making  use  of  (5.6),  (7.3)  and  (7.4),  we  find 


P  hi'  1  '  >  x 

(p  +  K(h))  A(h)R-(h)  <  lim  inf —L— — -1  (7.5) 

J  P  [aB{  >  x\ 

and  the  conclusion  (7.1)  is  an  immediate  consequence  of  (7.5)  and  of  the  equivalence  (6.6)  noted 
earlier.  ■ 


7.2  An  improved  lower  bound  for  two  sources 

This  time,  in  the  same  way  that  the  upper  bound  in  Proposition  6.1  leads  to  Proposition  6.2,  we 
have  the  following  result  from  Proposition  7.1. 


Proposition  7.2  Consider  a  superadditive  mapping  h  :  It+  —  It  such  that  whenever  0  <  e  <  e* 
for  some  e*  >  0,  the  scaled  mapping  eh  satisfies  (H2)-(H7)  [with  7  <—  c  —  p2  and  p  <—  pi].  Then,  it 
holds  that 


( lim  A  (eh) 
VT  0 


|limi?_(r/() 
V  v  1 0 


<  lim  inf 

x—>oo 


p  wAl+A2’c  >  x 
p  [WA1’C~P1  >  x] 


(7.6) 


The  lower  bound  (7.6)  can  be  improved  by  seeking  a  value  for  the  product  of  the  quantities 
liimjoA (eh)  and  liim jo  R-(eh)  that  is  as  large  as  possible  (if  not  the  largest)  among  admissible 
perturbations  h.  While  A  (eh)  <  1  is  always  true,  it  is  often  possible  to  argue  that 

lim  A(.e/i)  =  1  (7.7) 

40 

for  the  selected  h,  in  which  case  (7.6)  reads 


lim  R_(eh)  <  lim  inf 

£  J  O  x 


P  ir  '  • 

P  [II' '  r  P-‘  >  x] 


(7.8) 


For  instance,  under  the  condition  (5.9)  on  source  A  2,  we  find  A  (eh)  =  1  for  each  e  >  0  whenever 

(6.10)  holds  with  L  <  0.  If,  instead,  the  condition  (6.3)  holds  (L  =  0),  then  the  CLT  refinement 

(5.10)  to  (5.9)  needs  to  be  brought  into  the  picture.  If  h  is  a  convex  mapping  with  h( 0)  =  0,  then 
(6.3)  is  equivalent  to  h'(oo)  =  0,  and  the  requirements  (H2)-(H4)  naturally  lead  to  taking  h  =  —u 
for  some  concave  increasing  mapping  u  :  IR.+  — » IR.+  with  u'( 00)  =  0.  In  that  case  (5.11)  yields 


P[d2(t) 


p2t  >  eh(t )]  =  P 


Vt_ 

u(t ) 


z2(t)  >  -e 


t  >  0 


and  the  conclusion  A  (eh)  =  1  follows  from  (5.10)  if  —h  satisfies  (6.11)  with  H  =  00. 
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We  now  focus  on  finding  a  perturbation  function  h  such  that  lim.  n /?  '("//)  is  as  large  as  possible. 
With  h  =  —u  as  above, 


R-(sh)  =  lirninf 


P  [Bt  >  x] 


P  >  .r  —  o  1i?x(.r)] 


<  1,  <r>0, 


and  the  requirement  that  limfj0  R-(eh)  is  as  large  as  possible  will  be  met  if  lim£j oR-(sh)  =  1,  in 
which  case  the  bound  (7.8)  becomes 


1  <  lirninf 

x^oo 


p  WA1+AV  > 

P  [WAl'c~P2  >  x] 


Here  too,  whenever  B\  G  MDA( A),  Lemma  3.6  yields 


(7.9) 


R-(-£mBl) 


lim 

X—XXD 


P  [Bi  >  x ] 


P[H| 


>  x  —  a  lsmBl 


(*)] 


=  e 


f  >  0 


and  taking  h  =  — mBl ,  we  get  linn  j0  R-(sh)  =  1.  Thus,  in  most  cases  of  interest,  —mBl  is  expected 
to  be  the  perturbation  function  of  choice  for  getting  (7.9). 


8  Applications 


In  the  examples  which  we  now  discuss,  we  assume  that  the  fluid  A  is  the  superposition  of  two 
independent  fluid  processes  A\  and  A-j1  with  A\  a  stationary  “independent  on-off".  The  conditions 

Pi  +  P-2  <  c  <  ?’i  +  p-2  (8.1) 

are  enforced  throughout.  Propositions  6.2  and  7.2  will  be  invoked  in  that  context. 

8.1  In  the  maximum  domain  of  attraction  of  the  Gumbel  distribution 

As  indicated  already  in  some  of  the  comments  following  Propositions  6.2  and  7.2,  a  fairly  compre¬ 
hensive  discussion  should  be  expected  when  B\  G  MDA( A).  Before  indicating  in  Proposition  8.1 
the  extent  to  which  this  is  indeed  the  case,  we  present  the  following  technical  fact. 

Lemma  8.1  Consider  a  mapping  u  :  ]R.+  —  R+  which  is  strictly  increasing  and  concave  in  the 
limit,  i.e.,  there  exists  xo  >  0  such  that  the  restriction  u  :  [xo.oo)  — ►  It+  is  strictly  increasing  and 
concave.  There  exist  an  increasing  concave  mapping  U  :  It+  — >  It+,  and  a  constant  x *  >  xo  such 
that  U{ 0)  =  0,  Lr/(0+)  :=  lima,|0  U'(x)  is  finite  and  u(x)  =  U(x)  for  all  x  >  x*. 

Proof.  Define  the  mapping  u  :  R+  R+  by  u(x)  =  u(x o)  if  0  <  x  <  Xo  and  by  u(x)  =  u(x) 
if  x  >  xq.  This  mapping,  while  non-decreasing,  is  not  necessarily  concave.  Let  uc  denote  the 
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concave  hull  of  it,  i.e. ,  the  smallest  concave  mapping  uc  :  ]R+  — ►  ]R+  such  that  it,  <  iic.  In  fact, 
uc  =  —  conv(— it)  where  conv  denotes  the  convex  hull  operation  [26,  p.  36].  It  is  easy  to  check  that 
iic  is  increasing  and  concave  with  iic(x)  =  u(x)  from  some  x  onward.  The  desired  mapping  U  is 
now  obtained  by  taking  U{x)  =  min(:r,  uc[x))  for  all  x  >  0.  I 

The  function  U  associated  with  u,  through  Lemma  8.1  is  clearly  not  unique.  In  specific  examples  the 
one  constructed  in  the  proof  can  be  safely  replaced  by  a  more  natural  one  which  derives  naturally 
from  the  form  of  u . 

Proposition  8.1  Assume  B\  and  I\  to  be  1R ^-valued  rvs  with  0  <  E  [I?i]  ,  E[Ii]  <  oo,  such  that 
B\  G  C  and  Bf  E  S.  Suppose  that  the  function  mg1  given  by  (3.8)  is  asymptotically  equivalent 
to  some  mapping  u  :  1R.+  — *■  1R.+  which  is  strictly  increasing  and  concave  in  the  limit ,  and  let 
U  :  1R+  — *■  1R.+  denote  any  increasing  concave  mapping  associated  with  u  as  in  Lemma  8.1. 

If  B\  G  MDA(A),  then  the  following  holds: 

(1)  There  exists  a  subadditive  mapping  li  :  1R+  — *■  1R  such  that  whenever  0  <  e  <  e*  for  some 

s*  >  0,  the  scaled  mapping  eh  satisfies  (H1)-(H7)  [with  7  c  —  P2  and  p  pi],  A  possible 

choice  is  h  =  U ,  in  which  case 

lim  R+{sh)  =  lime"  =  1;  (8.2) 

fj.0  cJO 

(2)  There  exists  a  superadditive  mapping  h  :  1R.+  It  such  that  whenever  0  <  t  <  e*  for  some 
e*  >  0,  the  scaled  mapping  eh  satisfies  (H2)-(H7)  [with  7  <—  0—  p2  and  p  p\[.  A  possible 
choice  is  h  =  —U,  in  which  case 

lim R_(eh)  =  lime-"  lf  =  1.  (8.3) 

si 0  5 JO 

A  proof  of  this  result  is  available  in  Appendix  E. 

We  next  indicate  how  Proposition  8.1  applies  in  some  specific  cases.  It  is  worth  pointing  out  that 
in  the  context  of  Proposition  8.1,  the  task  of  checking  (6.7)  over  some  entire  interval  (0,e*|  can  be 

greatly  reduced:  Indeed,  with  h  =  U,  owing  to  the  tail  equivalence  of  af{B()  and  oB(  (0  <  e  <  e*) 

noted  in  the  proof  of  Proposition  8.1,  we  get  that  (6.7)  is  equivalent  to 

P  |  ir '  h  >  ,r]  =  O  (P  [aB\  >  x]) ,  0  <  e  <  £*.  (8.4) 

8.2  Log-normal 

By  definition,  we  have  B\  ,1.  ez  where  Z  is  normally  distributed  with  mean  //  and  variance  62 . 
Recall  that  B\  G  *S  (and  therefore  B\  G  C)  and  B\  G  S  [14,  Example  1.4.7,  p.  55].  It  is  also  well 
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known  that  B \  G  MDA(K)  [14,  p.  150].  Then,  straightforward  calculations  yield  the  asymptotic 
equivalences 

fi(  ;log:r  /,;7(2<G’) 


and 


Consequently, 


P[5l>'r]  "  \Z2tt  (log  x  —  //)  ’ 

,C!  .  ,  (l°g.r  Ilf  i 

P|B?>^E[Bl]^„og,-,,V  ^  “I*!  =“■»/+'•>• 

mLN(x)  ~  <52 


(8.5) 


(8.6) 


log  ;r  —  n  log  ;r 

with  m^jv  the  mapping  associated  through  (3.8)  with  the  rv  A'  =  B\.  Note  that  iuln  is  asymptot¬ 
ically  equivalent  to  the  mapping  uln  :  ]R.+  — ►  1R.+  given  by  uln(x)  =  t2{x/\og  x  V  e )  (x  >  0).  The 
mapping  uln  is  concave  and  strictly  increasing  in  the  limit,  whence  Proposition  8.1  can  be  invoked 
with  Uln  denoting  the  mapping  associated  with  uln  through  Lemma  8.1. 


Lower  bound:  Applying  the  lower  bound  (7.6)  with  h  =  —Uln  and  using  (8.3),  we  find 

P  [tr  '  • 


>  x 

limA(tL)  <  lirninf  „  rTT,  , 

•  :0  x-*va  P  [IP  '  '  P'2  >  x] 


(8.7) 


with 


A  (eh)  =  lim  inf  P 

t—>o O 


A2{t)  >  p2t  -  eb 


2  t 


s  >  0. 


log  t  J  ’ 

Assuming  the  CLT  refinement  (5.10)  to  hold  for  source  A 2,  we  conclude  that 

c2  v7! 


A  (eh)  =  lirninf  P 

t—>o O 


Z2(t)  >  —ebz 


log  t 


=  lim  P 


t~* oo  [ 


log  t 


Z2(t)  >  -eb2 


=  1, 


and  we  achieve  the  best  possible  lower  bound 


1  <  lirninf 


W 


Ai~i~A2,c 


>  X 


oo  P  [WA !’c  p2  >  x] 


(8.8) 


Upper  bound:  Applying  the  upper  bound  (6.8)  with  h  =  Uln  and  making  use  of  (8.2),  we  obtain 


lim  sup 

a?— >-oo 


p  ir  !  ■  >  x 

p  [ir  '  '  I"2  >  ,r] 


<  1 


(8.9) 


provided  the  second  source  A2  satisfies  the  condition  (6.7),  or  equivalently  (8.4).  We  now  dis¬ 
cuss  this  condition  when  A2  is  a  regenerative  source  as  defined  in  Section  6.3  under  the  moment 
conditions  (6.13).  With  the  choice  h  =  Uln ,  for  each  5  >  0  we  have  linq_.:x;:  eh(t)/l  =  0  and 


26 


Hindoo  sh{i)/tll2+v°  =  oo  for  each  v$  in  (0,1/2).  Therefore,  Proposition  6.3  applies  with  L  =  0 
and  H  =  oo,  we  are  lead  to  the  conclusion  that  for  each  t  >  0,  there  exist  finite  constants  0  <  6\,e 
and  6-2,s  >  0  such  that 


p  \yM,P2-sh  y  x 

limsup —  6  — -  <  82, s,  v  <E  (0, 1).  (8.10) 

x — *oo  xL~v  e  0l’sX 

It  is  now  plain  from  (8.6)  and  (8.10)  that  (8.4)  indeed  holds.  We  summarize  the  discussion  as 
follows: 

Proposition  8.2  Let  A  =  A\  +  A-j  be  the  superposition  of  two  independent  fluid  processes  A\  and 
A 2.  Assume  A\  to  be  a  stationary  “ independent  on-off”  source  such  that  pi  +  p-2  <  c  <  r\  +  p2  with 
activity  period  B\  distributed  according  to  (8.5),  and  let  A2  be  a  regenerative  source  (in  the  sense 
discussed  in  Section  6.3)  under  the  moment  conditions  (6.13).  Then,  it  holds  that 

P  111'  >  x 

1  i  121  - - - 

x  'Oo  P  (WAlX-p-i  >  x] 


The  lower  bound  (8.8)  only  requires  that  source  2  satisfies  (5.10).  Proposition  8.2  also  holds  when 
B\  has  a  Benktander-type-I  distribution  [14,  p.  149]  since  in  that  case  B\  G  MDA(A)  fi5  with  the 
corresponding  function  (3.8)  being  asymptotically  equivalent  to  kiln- 


8.3  Weibull 


With  a  >  0  and  0  <  0  <  1, 


P  [Hi  >  x]  =  exp  1^'  ax  ,  x  >  0. 


(8.12) 


Again  we  have  that  B\  e  S  fl  MDA(  A)  (and  therefore  B\  G  C  since  S  C  C)  [14,  Example  1.4.7,  p. 
55]  and  B*  G  5  [14,  p.  150].  All  the  moments  of  this  distribution  are  finite,  and  in  particular,  we 
have 

nr1) 


E  [B0  = 


0  a  p 


where  the  Gamma  function  P  :  (0,  00)  —  1R.+  is  defined  by 

POO 

T(s)  :=  /  xs~1e~xdx ,  s  >  0. 

Jo 

By  appealing  to  properties  of  the  incomplete  Gamma  function,  we  readily  see  that 


P  [Bf  >  x] 


1-/3 

a  13 

W 


m 


x1  J  exp 


i1'  ax  '  \  , 


(8.13) 
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and  we  find 


mw(x)  ~ 


1 


1-/3 


where  m\v  denotes  the  mapping  associated  through  (3.8)  with  the  Weibull  rv  A'  =  B \ 


Note  that  niw  is  asymptotically  equivalent  to  uw  ■  IR  -  — *■  1R+  given  by  uw(%)  '■=  ( a/3)~1x1~l } 
( x  >  0).  The  mapping  uw  is  concave  and  strictly  increasing  on  ]R+  with  Mm(0)  =  0,  and  Proposition 
8.1  can  be  applied  with  Uw  denoting  the  mapping  associated  with  uw  through  Lemma  8.1. 


Lower  bound:  Applying  the  lower  bound  (7.6)  with  h  =  —Lhv  and  using  (8.3),  we  also  find  (8.7) 
with 

£ 


A  (sh)  =  liminfP 


t—>o O 


8  >  0. 


Under  the  CLT  refinement  (5.10)  for  source  A 2,  we  see  that 


A  (sh)  =  liminfP 

t — >00 


zM  >  ~tl-> 


=  lim  P 

t — >00 


e-izM  > 


=  1 


provided  the  condition 


0  <  p  <  - 


(8.14) 


holds,  in  which  case  the  best  possible  lower  bound  (8.8)  is  achieved. 


Upper  bound:  Applying  the  upper  bound  (6.8)  with  h  =  Uw  and  making  use  of  (8.2),  we  obtain 
the  upper  bound  (8.9)  provided  the  second  source  A 2  satisfies  the  condition  (6.7).  This  condition  is 
discussed  now  when  A 2  is  a  regenerative  source  under  the  moment  conditions  (6.13).  For  the  choice 
h  =  Uw,  for  each  £  >  0  we  have  lim^oc  eh{t)/t  =  0  and  linp . .  x;  eh(t) /t1l2+v°  =  00  for  each  vq  in 
(0, 1/2)  such  that 

0  +  vO<\-  (8.15) 

Hence,  whenever  we  pick  v 0  in  (0, 1/2)  such  that  (8.15)  (thus  (8.14))  holds,  we  can  invoke  Proposition 
6.3  while  still  guaranteeing  (8.8).  Hence,  with  each  £  >  0,  there  exist  finite  constants  0  <  8%g  and 
62j£  >  0  such  that  (8.10)  still  holds  but  only  for  v  in  the  interval  (0, 1  —  20).  It  is  now  plain  from 
(8.13)  that  (8.4)  indeed  holds  if  there  exists  such  an  admissible  value  of  v  with  the  property  that 


lim  sup 

X^OO 


.r1  v  exp  (— bi^£xu) 
.r1--5  exp  (—aa~/5xd) 


=  0, 


(8.16) 


a  requirement  equivalent  to  0  <  v.  Therefore,  (8.4)  will  hold  if  there  exists  v  in  (0, 1)  such  that 
0  <  v  <  1  —  20 ,  a  non- vacuous  condition  only  if 

0  <  ^  (8-17) 


We  summarize  the  findings  as  follows: 
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Proposition  8.3  Let  A  =  A\  +  A2  be  the  superposition  of  two  fluid  processes  A\  and  A2 .  Assume 
A\  to  be  a  stationary  ‘independent  on-off”  source  such  that  pi  +  P2  <  c  <  r\  +  p2  with  activity 
period  B\  distributed  according  to  (8.12)  with  (8.17),  and  let  A2  be  a  regenerative  source  under  the 
moment  conditions  (6.13).  Then,  (8.11)  holds. 


Here,  the  lower  bound  (8.8)  requires  only  that  source  2  satisfies  (5.10)  but  under  the  additional 
condition  (8.14)  which  defines  the  so-called  moderately  heavy-tail  case.  On  the  other  hand,  the 
upper  bound  (8.9)  is  shown  to  hold  only  under  the  more  stringent  condition  (8.17);  this  constraint 
amounts  to  period  of  activity  B\  of  source  1  to  being  heavy-tailed  enough!  At  this  point,  the  reader 
may  wonder  whether  Proposition  8.3  still  holds  when  the  parameter  j3  in  (8.12)  lies  in  the  interval 
(1/3, 1).  Of  course,  such  a  conclusion,  if  correct,  would  have  to  be  reached  by  arguments  different 
from  the  ones  used  here.  The  following  fact  due  to  Dumas  and  Simonian  [11]  implies  a  partial 
negative  answer  to  the  question: 


Proposition  8.4  Let  A  =  A\  +  A2  be  the  superposition  of  two  independent  fluid  processes  A\  and 
A2.  Assume  A\  to  be  a  stationary  “independent  on-off”  source  such  that  pi  +  P2  <  c  <  r\  +  P2  with 
activity  period  B\,  and  let  A2  be  a  fluid  source  which  satisfies  (5.10).  Whenever 


lim  =  0 


(8.18) 


with  mBl  given  in  (3.8),  then  it  holds  that 


lim 

x—>oo 


P  \wA'c  >  x 

P  [VJMi,C-P2  > 


=  00. 


(8.19) 


Dumas  and  Simonian  establish  this  negative  fact  (8.19)  through  a  very  simple  argument  akin  to 
the  one  used  in  deriving  the  generic  lower  bound  in  Proposition  2.1.  When  applied  to  the  setup  of 
Proposition  8.3  with  ^  <  /3  <  1,  Proposition  8.4  implies  the  failure  of  the  equivalence  (8.11)  since 
we  now  have 


lim 

x—>oo 


mBl(x) 


lim  - — —  =  0. 

:r— 00 


A  condition  similar  to  (8.18)  was  also  encountered  in  recent  work  by  Asmusssen,  Kluppelberg  and 
Sigman  [3,  Thm.  4.1]  on  distributional  properties  of  the  sample  of  a  process  at  subexponential  times; 
there  as  well  connections  with  Extreme  Value  Theory  naturally  emerge.  When  specialized  to  the 
family  of  Weibull  distributions,  their  results  do  hold  for  (3  in  the  entire  range  (0, 1/2).  This  analogy 
holds  up  the  possibility  that  Proposition  8.3  might  indeed  be  valid  for  (3  in  the  interval  (1/3, 1/2). 
After  all,  the  argument  behind  the  upper  bound  in  Proposition  8.3  relies  in  an  essential  manner  on 
the  decay  rates  given  in  Proposition  6.3;  there  is  no  reason  a  priori  to  believe  that  they  are  best 
and  cannot  be  improved!  As  a  case  in  point,  we  remark  that  the  reduced  load  approximation  (8.11) 
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does  hold  when  the  rvs  { ( V',> .  Cn),  n  =  2, . . .}  characterizing  source  2  are  deterministic  (as  would  be 
the  case  for  an  on-off  source  .42  with  deterministic  on  and  off  periods).  Indeed,  we  then  get  (6.15) 

with  no  =  1/2  as  we  note  that  E  jV(V2_P2<:'2)]  =  1  in  (D.13)  so  that  g(t)2  jt  may  be  replaced  by  g(t) 
in  (D.18),  and  (8.16)  holds  with  v  >  (3. 

When  B\  has  a  Benktander-type-II  distribution  [14,  p.  149],  it  is  also  the  case  that  B\  E  MDA(  A)fl 
«F,  and  under  the  condition  (8.14),  Proposition  8.3  also  holds  since  the  corresponding  function  (3.8) 
is  asymptotically  equivalent  to  raw- 


8.4  Generalized  Pareto 


With  1  <  /3  <  2,  this  corresponds  to 

P  [ill  >■#].'=  x  ' /. ( .r ),  x  >  0  (8.20) 

for  some  slowly  varying  function  L  :  1R+  — ►  ]R.+  .  The  rv  B\  is  integrable  but  with  infinite  variance, 
and  its  integrated  tail  distribution  is  given  by 

P  [B\  >  x]  ~  E  B{  1  x  ;  1  /.(.r). 

We  denote  by  nip  the  mapping  associated  through  (3.8)  with  the  rv  A'  =  B\. 

Generalized  Pareto  rvs  (8.20)  do  not  belong  to  MDA(A)  but  to  the  maximum  domain  of  attraction 
of  the  Frechet  distribution  [14,  p.  121].  The  reader  will  find  this  matter  discussed  in  [14,  Thin. 
3.3.7,  p.  131].  Consequently,  Proposition  8.1  cannot  be  invoked,  and  while  this  may  be  viewed  as 
an  unfortunate  development,  we  shall  see  shortly  that  the  case  of  generalized  Pareto  rvs  is  in  fact 
easier  than  the  cases  treated  thus  far. 


Indeed,  we  have 


mp(x)  ~  x, 


and  in  the  spirit  of  Proposition  8.1,  nip  is  asymptotically  equivalent  to  the  (strictly  increasing) 
linear  mapping  up  :  x  — ►  x  (so  that  here  we  can  take  Up  to  coincide  with  up  in  Lemma  8.1).  This 
suggests  taking  perturbation  mappings  which  are  linear.  For  each  g  in  II,  the  mapping  gUp  does 
satisfy  conditions  (H1)-(H4)  provided 


c  -  (/'i  +  P2)  <  V  <  c  -  (pi  +  i>>). 

a  non-vacuous  constraint  under  (8.1).  With  the  notation  (E.2)  and  (E.3),  the  integrability  of  the 
rvs  B\  and  Ii  implies  that  of  the  rvs  A'/  and  av(B\)  (thus  of  a+(B i)).  Under  the  stability  condition 
we  then  find  E  [A’/]  <  0  for  g  in  a  small  enough  neighborhood  of  the  origin,  whence  (H5)  holds  for 
gUp.  Finally,  Bi  E  C  and  Bf  E  S  [14],  so  that  a+(B\)  E  C  and  a^(B*)  E  S  by  linearity,  hence 
(H6)  and  (H7)  are  satisfied.  Consequently,  for  g  in  a  small  enough  neighborhood  of  the  origin,  the 
mapping  gUp  does  satisfy  conditions  (H1)-(H7)! 
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Lower  bound:  We  apply  Proposition  6.2  with  h 
as  well  upon  noting 


A  (eh)  =  P 


'Mt) 

t 


=  —Up.  We  readily  conclude  that  (8.8)  holds  here 


p2  >  -£ 


1,  t  >  0 


under  the  condition  (5.9),  and  that 


R-{sh)  =  lirninf 


P  [(«  —  e)B\  >  .r] 
P  [ aB f  >  x] 


e  >  0. 


Upper  bound:  This  time  Proposition  6.2  will  be  applied  with  h  =  Up.  We  have 


R+(sh)  =  lim  sup 


P  [(a  +  s)Bf  >  x]  /n  ■  '  1 

P  (1  /i(  >  x]  V  OL 


s  >  0 


(8.21) 


so  that  (8.9)  is  also  obtained  provided  condition  (6.7)  holds.  Here  as  well,  this  condition  takes  the 
simplified  form 

P  |ir  >  x\  =  o  (P  [aB\  >  x])  ,  0  <  e  <  t*  (8.22) 

owing  to  the  fact  that  the  rvs  {(a  +  s)Bf,  e  >  0}  are  all  tail  equivalent  to  aBf  by  virtue  of  (8.21). 
If  the  second  source  A 2  is  a  regenerative  source  under  the  moment  conditions  (6.13),  then  (8.22) 
always  holds.  In  the  Pareto  case,  the  validity  of  (6.7),  or  equivalently  of  (8.22),  holds  more  widely, 
even  when  source  A 2  fails  to  have  finite  exponential  moments.  For  instance,  if  source  A 2  is  an 
independent  on-off  source  with  generic  activity  period  B-j  such  that  B?  G  C  and  B%  E  <S,  then  by 
Remark  4.1  we  have 

P  [wM,P2+s  >  a]  ~  k£ p  [(r2  -  (P2  +  s))BZ  >  .r] 

for  some  appropriate  constant  Ks  >  0  determined  by  the  source  statistics.  Hence,  (8.22)  holds 
provided  the  condition 

P  [(»’2  -  (P2  +  e))B2  >  x]  =  o(P  [u/l|‘  >  m] ) ,  0  <  e  <  e* 

is  met.  Collecting  all  these  remarks  leads  to  the  following  result  first  obtained  by  Jelenkovic  and 
Lazar  [18]. 


Proposition  8.5  Let  A  =  A\  +  A2  be  the  superposition  of  two  independent  fluid  processes  A\  and 
A'2 .  Assume  A\  to  be  a  stationary  “ independent  on-off”  source  such  that  pi  +  P2  <  c  <  r\  +  p2  with 
activity  period  B\  distributed  according  to  (8.20)  with  (5.9),  and  let  A2  be  a  fluid  source  satisfying 
(8.22)  (e.g.,  source  2  is  a  regenerative  source  under  the  moment  conditions  (6.13)).  Then,  (8.11) 
holds. 


We  complete  the  discussion  by  noting  that  if  A'  is  a  generalized  Pareto  rv,  then  for  any  mapping 
tp  :  IR.+  — >  IR.+  which  is  asymptotically  equivalent  to  mp,  we  have 


lim  p  [X  >  x  +  jMx)] 

5-+00  p  [A'  >  x] 


(1  +  y)  13  *  y  >  -1 


The  reader  will  note  the  analogy  with  a  similar  limit  in  Lemma  3.6  for  rvs  in  MDA(K ). 
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9  Additional  results  for  intermediate  regular  varying  rvs 


In  this  section  we  strengthen  Proposition  8.5  along  several  directions,  thereby  confirming  the  singular 
position  occupied  by  generalized  Pareto  rvs  and  their  extensions. 


9.1  Intermediate  Regular  Variation 


Jelenkovic  and  Lazar  [18,  Thm.  10]  show  that  Proposition  8.5  remains  true  if  B\  belongs  to  the 
larger  class  of  intermediate  regular  varying  rvs:  Following  [10,  Definition  (1.2)],  we  say  that  an 
R+-valued  rv  A'  is  an  interm, ediate  regular  varying  rv,  denoted  A  G  IB,  if 


lim  lim  inf 

AJ.1 


P  [A  >  Xx] 
P  [A  >  ,r] 


lim  lim  sup 

A  jl  x — ^oo 


P  [A  >  Xx] 
P  [A  >  x] 


=  1. 


(9.1) 


The  defining  relation  (9.1)  is  easily  seen  to  be  equivalent  to 


lim  lim  inf 

A 1 1 


P  [A  >  Xx] 
P  [A  >  ,r] 


lim  lim  sup 

A 1 1  x — too 


P  [A  >  Xx] 
P  [A  >  ,r] 


=  1. 


(9.2) 


It  is  known  that  B  C  IB  C  S,  where  B  is  the  class  of  regular  varying  distributions  (which  coincides 
with  the  class  of  generalized  Pareto  rvs)  [5,  p.  18]. 


Our  approach  can  also  be  used  to  extend  the  validity  of  Proposition  8.5  to  the  case  when  B\  G  TR. 
The  proof  is  analogous  to  that  of  Proposition  8.5  upon  selecting  the  mapping  h(x)  =  —x  in  the  lower 
bound  and  h{x)  =  x  in  the  upper  bound.  The  only  difference  is  that  now  the  limits  lime jq  R-  (f-h )  =  1 
(with  h(x)  =  —x)  and  limejo  R+(eh)  =  1  (with  h(x)  =  x)  follow  directly  from  (9.1)-(9.2)  together 
with  the  (easily  checked)  property  that  B*  G  IB  if  B\  G  IB  with  0  <  E  [B\]  <  oo. 


9.2  Upper  bound  without  independence 

This  last  result  obtained  by  Jelenkovic  and  Lazar  [18,  Thm  10]  can  be  improved  along  yet  another 
direction. 


Proposition  9.1  Let  A  =  A\  +  A 2  be  the  superposition  of  two  fluid  processes  A\  and  A^.  Assume 
A\  to  be  a  stationary  “independent  on-off”  source  such  that  pi  +  P2  <  c  <  r\+p2  with  activity  period 
B\  G  IB  and  let  A2  be  a  fluid  source  satisfying  (8.22)  for  some  s*  >  0.  Then,  it  holds  that 


lim  sup 

X^OO 


p  fir  '  r  >  .<•' 

P  [U/V4i,c-p2  >  _r] 


<  1. 


(9.3) 


In  contrast  with  Theorem  10  in  [18]  we  do  not  require  sources  A\  and  A2  to  be  independent.  This 
is  made  possible  by  the  asymptotic  scale  invariance  implied  by  (9.1)  and  (9.2). 
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Proof.  Pick  £  in  (0,  min(c  —  (pi  +  P2),-*))  as  in  (8.22).  Lemma  2.1  with  h{t )  =  st  (t  >  0)  and  a 
standard  union  bound  argument  yield 

P  [if  '•?  >  a;]  <  P  \wAl'c~p2~s  >  ?/x]  +  P  [wM’P2+£  >  (1  -  r])x\  ,  x  >  0  (9.4) 

for  any  rj  in  (0, 1).  Using  (4.12)  [with  7  <—  c  —  p-2  —  s,  p  <—  pi]  for  B\  G  Z1Z  (thus  B\  G  HZ),  we  get 

p  '  •  >  rjx] - - p  +  £)_g+  >  > .  v  g  [0i  1])  £-  G  [0,  £*)  (9.5) 

L  J  c  ~(pi+  P2)  ~  e 

with  a  given  by  (5.1)  as  usual. 


By  writing 


p 

WA 2,P2+£  >  (1  _  j-j'jx 

P 

ir  1  ■  >  (i  —  i])x 

p 

WA i,c  p 2  >  (I  //)./■ 

P  [WAl’c~p2  >  x] 

"  ““  ~P 

\wAuc-P2  >  (i  -  n)x] 

P  [J,p.4i,c-p2  >  x] 

we  conclude  from  (8.22),  (9.2)  and  (9.5)  that 


lim  lim  sup 

r)]  1  x^oo 


p  [ir  xi 

P  [  1U J  b  ; C  P 2  >  .p] 


=  0. 


On  the  other  hand,  application  of  (9.5)  on  each  of  the  factors  in 


p  U"  L  '  P2  £  >  rjx  P  WAl,c  P2  s  >  r)x  P  U"  '  '  p 2  >  7)x 

p  [p/-b,c-p 2  >  p  [ppXc-pa  >  jp]  p  [VU^1,C-P2  >  x] 


gives 


lim  sup 

X^OO 


p  (L  !  '  f2  E>  jjx 
p  ir  '  '  >  x] 


c~(pi  +  p2)  ..  P  [iff  >  a(a  +  g)  \r] 

C-(P1+P2)-S  ,  .x  1  P  [iff  >  a1] 


lim  sup 

x—>oo 


P  [iff  >  ?/;r] 
P  [iff  >  x]  ' 


Hence,  by  appealing  twice  to  (9.2)  and  using  (9.7)  we  obtain 


p  \yAi$  p2  s  >  7jX 


lim  lim  lim  sup  ■  „  , 

£J0  p|l  a--oc  P  [IT  '  '  p-  >  X 


<  1, 


and  combining  now  (9.4),  (9.6)  and  (9.8)  readily  leads  to  (9.3). 


(9.6) 


(9.7) 


(9.8) 


10  Superposition  of  independent  fluid  sources 


Sections  5  onward  have  dealt  with  the  multiplexing  of  two  independent  fluid  sources  A\  and  A2, 
where  A\  was  assumed  to  be  a  stationary  “independent  on-off”  source  and  A2  was  arbitrary.  In  this 
section  we  consider  the  case  when  A\  (respectively,  A2)  is  in  turn  obtained  by  the  superposition  of 
independent  on-off  (respectively,  fluid)  sources. 
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10.1  Superposition  of  independent  on-off  sources  in  Ai 


We  first  consider  the  situation  where  the  fluid  process  A\  =  {Ai(t),  t  >  0}  results  from  the  su¬ 
perposition  of  N  mutually  independent  fluid  processes  A\$  :=  {Wg^f),  t  >  0},  i  =  1,2,...,  iV,- 
namely 

N 

Ax(t)  :=J2  AiAt),  t>  0.  (10.1) 

i= 1 

For  each  i  =  1 , 2, N',  the  source  A\i  is  assumed  to  be  a  “stationary  independent  on-off"  source 
with  peak  rate  r\  and  we  set  r±  :=  J2iLi  ri i- 

It  is  natural  to  seek  an  extension  of  the  reduced  load  approximation  result  that  follows  from  Propo¬ 
sitions  6.2  and  7.2  when  A\  is  a  superposition  of  on-off  sources  as  defined  above  and  r\  >  c  —  p2- 
Unfortunately,  such  an  extension  turns  out  to  be  extremely  difficult,  and  we  contend  ourselves  with 
only  a  lower  bound.  Proposition  10.2  generalizes  a  result  due  to  Choudhury  and  Whitt  [8,  Thm  3] 
by  allowing  the  presence  of  a  “background”  source  A2  that  essentially  reduces  the  service  capacity 
c  by  its  mean  rate  p2-  We  present  an  intermediate  result  first. 

For  each  i  =  1, . . . ,  N,  we  denote  by  {(In$i  Bn$,  n  =  0.1,...}  the  alternating  sequence  of  off  and  on 
periods  for  the  stationary  version  of  the  independent  on-off  source  A\^  as  described  in  Section  4.1. 


Proposition  10.1  Let  A\  =  J2tL  1  A\^  be  as  in  (10.1),  and  assume  r\  >  c  —  p2-  Then,  for  any 
mapping  h  :  1R.+  — *■  HI ,  ,  with  a(t)  :=  (rq  —  (c  —  P2))t  —  h(t )  (t  >  0),  it  holds  that 


N  ( 

[’ wM'C-p2'-h  >  .r]  >  n 

L  ;  1  V 


E[5 


l,d 


y  V®  I-®!,*]  +  E  [I] 


l,d 


a(Bf  t)  >  x 


x  >  0. 


(10.2) 


Proof.  Define  B(  :=  min,:=i  ...  jv -Bo  i  where  :=  BqA  [Tq,*  =0],  i  =  1 , . . . ,  IV.  Then,  starting 
from  the  obvious  bound 

MB’0)  -  (c  -  p2)B'0  -  h(B'0)  < 
and  noting  that  A(B'0)  =  r\B'0,  we  see  that 


P  |  ir  '  /••••  h  >  ,r]  >  p  [riB'0  -  (c  -  p2)B'0  -  h(B'0)  >  x] 

=  P  [a(B'0)>x] 

>  P  |  min  alB'nf)  >  x 

i= 1 . N  u,t/ 


=  n p  [«(/o  > 

i=  1 
N 

>  n  P  =  °]  P  [a(B0,i)  >  %\Io,i  =  0] 


i= 1 
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and  (10.2)  follows  from  the  fact  that  [a(-Bo,i)|-fo,i  =  0]  =st  a(Btj),  i  =  !,  •  •  • ,  N.  | 

The  next  result  follows  from  Proposition  10.1  and  Proposition  2.1  in  a  manner  similar  to  the  proof 
of  Proposition  7.1. 


Proposition  10.2  Let  Ai  =  J2iLi  Ai ,i  be  as  in  (10.1),  and  assume  r i  >  c  —  p2.  Then,  for  any 
mapping  li  :  IR+  — *■  II which  satisfies  (H4)  [with  r  <—  r\,  7  <—  c  —  p2,  h  < - h],  it  holds  that 


lim  inf 

X — ^00 


p 

ir  1  -  >  x 

nili  p ; 

(ri  -  c  +  p2)B\ 

i  >x) 

N 


N 


>A(h)l[R-4h)l[ 


E  \B 


hd 


i= 1 


+£[/!,] 

with  A (h)  given  by  (5.8)  and  R-j(h )  by  (5.6)  with  B\  replaced  by 


When  the  on  periods  i  of  .1 1  ^  have  generalized  Pareto,  log-normal,  or  Weibull  (with  0  <  /3  <  1/2) 
distributions  and  A 2  satisfies  a  CLT,  we  can  obtain 

N 

A(h)'[[R-<i(h)  =  l 

i= 1 

as  done  in  Section  8  by  chosing  an  appropriate  sequence  of  h.  In  that  case  we  obtain  the  same 
constant  as  [8,  Theorem  3].  More  generally,  by  choosing  h  =  0,  and  by  assuming  that  A 2  satisfies 
a  CLT  we  can  obtain 

N 

A{h)  JJ  R-ti(h)  =  1/2 

i= 1 

with  no  further  assumptions  on  the  on  periods  B\j  of  A\  i. 

10.2  Superposition  of  independent  fluid  sources  in  A2 

Propositions  6.2  and  7.2  require  fairly  mild  conditions  on  source  A2.  In  practice,  it  is  often  the  case 
that  source  A2  is  the  superposition  of  independent  fluid  processes,  say  regenerative  fluid  sources,  or 
even  more  specifically,  independent  on-off  sources.  The  question  thus  naturally  arises  as  whether 
the  requisite  conditions  on  the  aggregate  source  A2  which  appear  in  Propositions  6.2  and  7.2  are 
implied  by  these  conditions  on  the  component  sources. 

We  investigate  these  issues  in  the  following  context:  Let  A2y  :=  { A2j(f).  t  >  0},  i  =  1,2,....,  N, 
be  N  mutually  independent  fluid  processes  with  average  rates  ^2,1,  •  •  • ,  P2,n,  respectively.  The  fluid 
process  resulting  from  the  superposition  of  .  I2.1  • . . . ,  A2, n  is  the  fluid  process  A2  :=  { A2(f).  t  >  0} 
defined  by 

N 

A2(t)  :=  ^2  A2 M,  t  >  0;  (10.3) 

1=1 
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its  average  rate  is  given  by  p2  ■=  Sill  P2,i ■  As  we  refer  to  the  examples  treated  in  Sections  8  and  9, 
we  note  that  the  choice  of  appropriate  perturbation  directions  h  in  the  lower  and  upper  bounds  are 
governed  by  the  distribution  of  the  activity  period  of  source  1,  with  the  selection  ensuring  the  largest 
and  smallest  possible  values  for  liimjo  R-(sh)  and  liimjj  i?+(r/(),  respectively,  among  admissible 
perturbations.  It  remains  therefore  to  explore  how  the  impact  of  the  individual  sources  on  these 
bounds  affects  the  impact  of  the  aggregate  source  A-2 . 


Lower  bound:  In  Proposition  7.2,  the  contribution  of  source  A2  to  the  lower  bound  (7.6)  arises  only 
through  the  constant  limfjo  A(sh)  for  some  appropriate  perturbation  function  h  :  IR.+  —  IR,  and 
it  is  desirable  to  have  (7.7),  i.e.,  limfjo  A(e/i)  =  1.  The  mutual  independence  of  the  component 
sources  ^2,1 , . . . ,  y  implies 


N 


A,s '{eh)  :=  TT  liminfP 


;=i 


1  -  P2,it  >  —h(l  ) 


<  A  (eh),  e  >  0 


so  that  we  may  substitute  the  constant  lining  As(eh)  in  the  left-hand  side  of  (7.6)  for  liimjo  A  (eh). 
Hence,  if  for  some  5*  >  0,  we  have 

Ai(sh)  :=  liminfP  [A^f)  —  p2pt  >  sh(t )]  =  1,  0  <  s  <  s*  (10-4) 

i'foo  ’  ’ 

for  all  7  =  1,2,...,  N ,  or  even  simply, 

A*  lim  Ai(eh)  =  1,  1  =  1,2 . N,  (10.5) 

fj.0 

then  limfjo  Ag(eh)  =  1,  thus  liin£io  A(s// )  =  1.  In  other  words,  the  desired  requirement  (7.7)  on 
A2  is  implied  by  the  similar  requirement  (10.5)  on  each  of  the  sources  A 24, . . . ,  A2,n-  As  pointed 
out  earlier,  and  as  further  discussed  in  Section  8,  (10.4)  or  (10.5)  will  hold  in  many  cases  of  interest 
when,  depending  on  /i,  either  the  Law  of  Large  Numbers  (5.9)  or  the  Central  Limit  Theorem  (5.10) 
holds  for  each  of  the  processes  A2J, . . 


Upper  bound:  We  now  turn  to  the  upper  bound  (6.8)  in  Proposition  6.2  when  A2  is  given  by 
the  superposition  (10.3).  In  that  case,  the  required  condition  (6.7)  with  respect  to  some  mapping 
h  :  IR  :  —  IR.  will  read  as 


P 


/  N 

sup  (Afoit) 

L*>°  V,  1 


P2,it) 


>  X 


=  o(P  [af(Bt)>x]), 


0  <  5  <  £*. 


(10.6) 


Therefore,  if  the  conditions 

P  h  >  x /AT]  =  o  (P  [«g" (Bi)  >  x] )  ,  0  <  e  <  s*, 

for  all  i  =  1,  2, . . . ,  iV,  are  simultaneously  satisfied,  then  a  standard  union  bound  argument  implies 
that 


P 


>  x 


N 

Y  I’  :■>'  >  X/N^ 

i= 1 

=  o(P[af(Bt)>x]) 
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for  0  <  e  <  t*,  and  (10.6)  holds.  This  argument  does  not  require  that  the  sources  be  independent. 


11  Conclusions  and  open  problems 


Although  we  have  succeeded  in  providing  some  conditions  under  which  the  reduced  load  equivalence 
(1.2)  holds,  the  picture  is  far  from  complete,  with  many  questions  still  left  unanswered.  We  review 
some  of  them  below: 

In  view  of  the  negative  result  of  Dumas  and  Simonian  [11]  described  in  Proposition  8.4,  the  equiva¬ 
lence  (1.2)  cannot  hold  under  (8.18),  and  it  is  natural  to  speculate  as  to  the  form  of  the  asymptotics 

for  the  tail  probabilities  P  [w4l+"42’c  >  xj  even  in  the  simple  case  when  A-2  is  an  exponential  on-off 
source. 

Underlying  the  discussion  presented  here  is  the  non-triviality  condition  ?q  +  p2  >  c  to  ensure  that 
source  1  is  not  immediately  flushed  out  when  fluid  is  released  at  the  reduced  rate  c  —  p2-  The 
equivalence  (1.2)  is  therefore  meaningless  if  r\  +  p-2  <  c,  and  a  completely  different  approach  is 

needed  for  obtaining  the  correct  asymptotics  of  P  [lb  '  '  >  xj  even  in  the  simple  case  when  A 2 

is  an  exponential  on-off  source. 

The  generic  condition  (6.1)  is  a  natural  one  for  establishing  the  upper  bounds  in  the  context  of 
subexponential  distributions,  e.g.,  Lemma  3.1(3)  and  the  line  of  argument  flowing  from  the  bound 
(6.4).  In  Section  6.3,  by  the  intermediary  of  Proposition  6.3,  we  are  now  in  possession  of  conditions 
to  check  the  validity  of  (6.1)  in  terms  of  the  rate  of  growth  for  the  perturbation  function  h  and 
the  statistics  of  the  source  A 2.  Proposition  6.3  (with  linear  perturbation  functions)  implies  that 
the  class  of  regenerative  on-off  sources  constitutes  a  subclass  of  the  class  of  exponential  sources 
introduced  by  Jelenkovic  and  Lazar  [18]  to  ensure  their  version  of  (6.1).  However,  what  was  needed 
here  is  an  estimate  on  the  rate  of  decay  of  tail  probabilities  associated  with  WA2'p2'~h  for  non- linear 
perturbation  functions!  In  establishing  this  rate  of  decay,  finite  exponential  moments  were  essential 
for  allowing  the  repeated  use  of  Chernoff  bounds.  Therefore,  several  questions  suggest  themselves 
very  naturaly:  As  the  discussion  following  Proposition  8.3  clearly  indicates,  a  better  decay  rate  is 
needed  if  one  is  to  handle  successfully  the  moderately  ligh  tailed  case  in  its  entirety.  Also,  it  is 
of  interest  to  find  out  what  happens  when  source  2,  while  still  regenerative,  does  not  have  finite 
exponential  moments;  a  completely  new  approach  would  be  required  to  get  the  appropriate  version 
of  Proposition  6.3  in  that  case. 

Finally,  in  the  introduction  we  mentionned  the  possibility  of  using  the  reduced  load  equivalence  (1.2) 
for  computational  purposes.  It  seems  intuitive  that  the  heavier  the  tail  of  B\  the  better  should  the 
approximation  be,  but  further  work  is  required  to  confirm  this  fact. 
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A  Proofs  of  Lemmas  3. 1-3. 5 


A  proof  of  Lemma  3.1. 


The  proof  of  Claim  1  is  straightforward  when  A'  G  C  and  available  in  [24,  Lemma  2]  when  A'  G  S. 
The  proof  of  Claim  2  is  an  easy  exercise  based  on  bounded  convergence  and  is  therefore  omitted. 


When  the  rv  T  is  non-negative,  the  proof  of  Claim  3  can  be  found  in  [4]  or  in  [9,  Thm.  1,  p.  533]. 
When  T  is  an  R- valued  rv,  the  proof  of  Claim  3  proceeds  as  follows:  Fix  x  >  0  and  note  that 

P  [A  +  T  >  x]  =  P  [A  +  T  >  x  \  T  >  0]  q  +  P  [A  —  |T|  >  x  \  T  <  0]  (1  -  q)  (A.l) 

with  q  :  =  P  [T  >  0].  Conditionally  on  [ T  >  0]  the  rvs  A'  and  T  are  independent  with 

P  [ T  >  x  |  T  >  0]  =  ^P  [T  >  x]  ^P  [Z  >  x] 

and 

P  A  >  x  /  >  0  ~  P  [A'  >  x]  ~  ciP  [Z  >  x]  . 

Thus,  applying  Lemma  3.1(3)  to  the  conditional  rvs  [A'|T  >  0](=sj:  X)  and  [T\T  >  0]  which  are 
non-negative,  we  find 

P  [A  +  T  >  x  |  T  >  0]  ~  (c,  +^P[Z>  .r]  . 

In  a  similar  way,  conditionally  on  [ T  <  0]  the  rvs  A'  and  T  are  independent  rvs  with 

P  A  >  x  ':  /  <  0  ~  P  [A'  >  x]  ~  ciP  [Z  >  x]  . 

Consequently,  by  Lemma  3.1(1),  the  rv  [A’  |  T  <  0]  belongs  to  5,  thus  to  £,  while  [|T|  |  T  <  0]  has 
support  in  R+.  A  straighforward  application  of  Lemma  3.1(2)  to  these  conditional  rvs  yields 


P  [A  -  |T|  >  x  |  T  <  0]  ~  P  [A  >  ,r]  ~  qP  [Z  >  x]  . 
The  proof  is  completed  upon  collecting  (A.1)-(A.2). 

A  proof  of  Lemma  3.2. 


It  is  plain  from  (3.2)  that 


P  [A  >  x] 
P  [A*  >  .r] 


E[A] 


[A'  >  x  +  t] 
[A  >  x] 


x  >  0. 


(A- 2) 


(A.3) 
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Because  A’  G  C.  we  also  have  lim, _ ^ 


P[-Y>a:4-tl 

P[X>x] 


1  for  each  t  >  0,  whence 


lim 

x—>oo 


P  [A'  >  x  +  f] 
P  [A  >  ^ 


dt  =  og 


(A-4) 


by  Fatou’s  lemma,  and  the  conclusion  (3.3)  is  now  a  straightforward  consequence  of  (A.3)-(A.4). 
Next,  for  each  x  >  0,  we  note  that 


E  [A]  P  [A*  >x]~  P  [A  -  Y  >  u]  du  =  /  (P  [A  >  u]  -  P  [A  >  u  +  Y]  )du 

J  X  J  X 

P  [A  >  u  +  Y] 


P  [A  >  u]  1  - 


P  [A  >  u] 


du. 


By  Lemma  3.1(2),  A'  and  A'  —  Y  have  the  same  right  tail,  i.e.,  for  every  c  >  0,  there  exists 
u*  =  u*{s)  >  0  such  that 


Consequently, 


p  a  >  u  +  y 

0  <  1 - i — -  <  S,  U  >  V  . 

P  A  >  u]  ~ 


1  r°°  P  [A  -  Y  >  u]  du 

0  <  1 - r-7~ - ““T"— - — ! -  <  £,  x>  u* 


E  [A]  P  [A*  >  x 
and  the  conclusion  (3.4)  immediately  follows. 


A  proof  of  Lemma  3.3. 


Under  the  assumptions  on  <p,  there  exists  ;r*  >  0  such  that  on  the  interval  [s*:,oo),  ip  is  strictly 
increasing  with  the  sets  </3([^*,co))  =  [y?(x*),oo)  and  (p([0,x*))  being  non-intersecting.  Moreover, 
x*  can  always  be  selected  large  enough  so  that  x*  >  xo  and  <p(x*)  >  0.  Consequently,  the  restriction 
of  p  to  [x*,  og)  is  a.e.  differentiable  and  invertible  with  liiiiy^^  p~1(y)  =  oo  and  we  have  {y  G  ]R.+  : 
■p(  //)  >  u}  =  (</3_1(w),  oo)  as  soon  as  <p(x*)  <  u. 

For  x  >  </?(x*),  we  note  that 


e  i  ^  •  ( -v  i ;  p  [(u+(A)r  >  *] 


P  [<yc(A')  >  u]  du 


P  [^?( A)  >  p{v)]  tp'{v)dv 


P  [A'  >  v]  p'(v)dv. 


(A- 5) 


The  lower  bound 


P[(^+(A))*>.r] 
P[A*  >  p-^x)] 


>  infjtp^x) 


v  >  p  1(x)} 


E[A] 

Eb+(A)] 
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is  immediate  and  letting  x  go  to  infinity  in  it,  we  conclude  that 


r  .  eP[&t(X))*>x]  ^  ,,  ,  ^  „  E[A] 

Inn  mi  _  r  — - — — A  >  Inn  ini{</3  (v)  :  v  >  tp  {x))-_ 


11111  1111  P  ,  v  i  11111 

x  P  [A*  >  tp~1(x) \  x->'ao 

Similarly,  we  have  the  upper  bound 

P[(^+(A))*>;r] 


Eb+(A)] 


so  that 


lim  sup 


P  [A*  >  tp-ifr)] 
P[i<p+(X))*>x] 


<  sup{  p'(  r)  :  v  >  tp  (./')} 


E[A] 


<  lim  sii|)^.',( r)  :  v  >  p  1(j)} 


E  [<P+(X)Y 

E[A] 


(A-6) 


P  [A'*  >  £>-!(*)]  '"A  "r  1  r  ^  ^  E  [</?+(  A)]  ’ 

The  result  (3.7)  readily  follows  from  (A. 6)  and  (A. 7)  under  the  existence  of  the  limit  (3.6). 

A  proof  of  Lemma  3.4. 


(A- 7) 


We  first  establish  the  result  in  the  special  case  when  the  mapping  tp  :  ]R.+  — ►  ]R.  is  strictly  increasing 
and  convex  on  ]R.+  with  <^(0)  =  0.  Under  these  assumptions,  tp  is  continuous  and  has  a  uniquely 
defined  inverse  tp -1  on  R+  with  </?_1(0)  =  0.  The  convexity  of  tp  implies  the  concavity  of  tp-1,  and 
we  have  lim.,. _ ^  tp(x)  =  liinx _ ^  <y2_1(:r)  =  oo. 


(Claim  1):  It  is  always  the  case  that 


lim  sup 

X^OO 


P  [tp(X)  >  x  +  y] 
P  [tp(X)  >  x] 


<  1, 


y  >  o. 


On  the  other  hand,  the  concavity  of  tp  1  coupled  with  tp  1(0)  =  0  implies  the  subadditivity  of  tp  1, 
i.e. ,  tp~1(x  +  y)  <  v?_1(.r)  +  tp~1(y)  for  all  x,  y  >  0.  Hence,  fixing  y  >  0,  we  find 


PMA  )>x  +  y]  =  PfAx^V  +  l/)] 

>  P[A  >p,-\x)  +  p-\y)\, 


x  >  0. 


Consequently, 

P  M  A)  >  ,r  +  y\  P[X  >p-1(x)+p>-1(y)] 
P  [p(A)  >  x]  ~  P  [A  > 

and  using  the  fact  that  A  G  C,  we  see  that 


P  [tp{  X )  >  x] 


>  1, 


y  >  o. 


The  desired  conclusion  tp( X)  G  C  follows. 
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(Claim  2):  This  time,  the  convexity  of  p  and  the  condition  <^(0)  =  0  imply  tp  superadditive,  i.e. , 
p(x  +  y)  >  p(x)  +  p(y)  for  all  x,  y  >  0.  Consequently,  with  a  rv  X'  distributed  like  X  but 
independent  of  it,  we  have 


so  that 


P  [p(X)  >  x 

The  condition  X  €  5  yields 


P  [ip(X)  +  p(X')  >  x]  <  P  [p( X  +  X')  >  x]  ,  x  >  0 


P  [<p(X.)  +  <p{X')  >  x]  <  P  [A"  +  X'  >  y-1(.r)]  ^  >  Q 


P  A  >  r-  '(.r) 


P  WA’)  +  p(X')  >  x\  _ 

Inn  sup - —7 — - - -  <  2, 

x^oo  P  [<p(A )  >  x] 

and  the  conclusion  p(X)  E  S  is  now  immediate  once  we  note  that  it  is  always  the  case  that 


x  >00  P  )  >  x\ 


>  2. 


We  now  turn  to  the  general  case  by  considering  a  mapping  p  which  satisfies  the  weaker  assumptions 
of  the  lemma.  By  convexity,  liin.,,^^  p(x)  =  00,  and  by  the  finiteness  of  <f>o,  there  exists  x *  >  Xq  such 
that  on  the  interval  [x*,oo),  tp  is  strictly  increasing  and  convex  with  <J?o  <  <p(x*)  and  0  <  ip(x*) 
(hence  0  <  tp[x)  for  x  >  x*).  Now,  consider  the  interpolated  mapping  ip*  :  1R.+  ]R.+  given 

by  <p*{x)  =  ^pAx  (0  <  x  <  x *)  and  <p*(x)  =  p{x)  (x*  <  x).  By  construction,  the  mapping 
p*  :  1R+  — >  R  is  strictly  increasing  and  convex  on  R+  with  y)*(0)  =  0.  Therefore,  by  the  first  part 
of  the  proof,  whenever  A'  belongs  to  C  (resp.  to  5),  it  follows  that  p*(X)  is  an  element  of  C  (resp. 
of  5).  The  desired  conclusion  on  r'(A )  now  follows  from  Lemma  3.1(1)  once  we  observe  that  the 
rvs  p(X)  and  p*(X)  have  ecjuivalent  right  tails,  i.e., 


lim 

X^OO 


P  b(A')  >  x] 
P  [p*{X)  >  x] 


1, 


a  fact  readily  verified  from  the  construction  of  h*  under  the  assumptions  on  h.  I 

A  proof  of  Lemma  3.5.  Fix  x  >  0  and  y  in  1R..  The  desired  conclusion  (3.10)  is  an  immediate 
consequence  of  the  obvious  relation 


P  [A'*  >  x  +  yp(x)\  m  \  (x  +  y p(x))  P  [A'  >  x  +  yp(x)\ 
P  [A'*  >  x]  ni\  (x)  P  [A'  >  x] 


and  of  the  assumed  limit  (3.11). 
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B  A  proof  of  Proposition  4.1 


We  drop  the  superscripts  A,  h  and  7  from  the  notation  as  these  quantities  are  fixed  throughout  the 
discussion.  As  in  Section  4.1  we  write  Tq  =  0  and  Tn+ 1  :=  Ylk=o (4  +  Bk)  for  n  =  0, 1,  . . ..  We  also 
note  that  A(Tn+ 1)  =  r  0  Bk  for  n  =  0, 1, . . .. 

(Claim  1):  Restricting  attention  to  the  epochs  {Tn,  n  =  0,1,...}  in  the  supremum  entering  the 
definition  of  W  and  noting  that 


n—  1 

h(Tn)  >  J2  h(h)  +  h(Bk), 


k= o 


n  =  1,2,... 


(B.l) 


by  the  superadditivity  of  h,  we  obtain  the  bound 

W  >  sup  {A{Tn)  —  ^Tn  +  h(Tn)) 

n= 0,1,... 


>  max  ^/i(0),  sup  (rRfc  -  y(4  +  4b)  +  M4)  +  h(Bk))J  j 


=  max  |  4(0) ,  sup  (  E  A'fc  =  V. 


ir.  -o.i.... 


Vfc= 0 


(Claim  2):  Condition  (HI)  (resp.  (H2))  ensures  that  the  mapping  a  (resp.  b)  is  monotone  increasing 
(resp.  decreasing),  and  it  is  therefore  easy  to  check  that 

W  =  sup  (A(Tn)  -  iTn  +  h(Tn)%  (B.2) 

n=0,l,... 

where  this  time  around,  the  subadditivity  of  h  yields 


n—  1 

h{Tn)  <  E  h(Ik)  +  h(Bk ),  n  =  1,  2, -  (B.3) 

k= 0 


Consequently,  by  an  argument  similar  to  that  given  in  the  proof  of  Claim  1,  we  get 


W  <  max  h{  0),  sup 

V 

=  max  /?,(0),  sup 

\  n= 0,1, 


^E  44  ~  7(4  +  Bk)  +  h(Ik)  +  /i(4))H 
4=o  ) ) 


~(£Avrv‘ 
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C  A  proof  of  Proposition  4.2 


Here  as  well,  we  drop  the  superscripts  A,  h  and  7  from  the  notation.  Note  from  (4.9)  the  relation 

P  [V  >x]=P  [X0  +  M  >  x] ,  x  >  h( 0)  (C.l) 

with  A'o  independent  of  the  rv  M  which  is  given  by 

M:=(  sup  .  (C.2) 

\n=lX...k=1  J 

As  we  have  in  mind  to  invoke  Lemma  3.1(3),  we  consider  in  turn  the  asymptotic  behavior  of  each 
the  rvs  A'o  and  M. 


(Step  1)  The  discussion  will  make  use  of  several  technical  facts  which  we  now  develop:  Under  (H2) 
the  inequality- 


holds,  so  that 


and 


h(t)  —  h{ 0)  =  [  h'(s)ds  <  7t,  t  >  0 

J  o 


6(6)  <  b( 0),  t  >  0 


(C.3) 


a+{t)  <  (( /•  —  y)t  +  h( 0)  +  7 1)+  <  h( 0)+  +  rf,  t  >  0.  (C.4) 

As  a  result  of  (C.4),  the  integrability  of  B\  implies  that  of  a+{Bi);  Tchebychev’s  inequality  now 
yields 

,rP  [a+(Hi)  >  x]  <  E  [a+{Bi)]  ,  x  >  0  (C.5) 

and  the  rv  a+{B i)  being  long-tailed  by  (H6),  we  can  select  x  large  enough  to  conclude 

0  <  E  [a+(Bi)]  <  oo.  (C.6) 

Next,  under  (H3)  and  (H4)  it  follows  from  (C.4)  that  a+  satisfies  the  conditions  (i)  and  (ii)  of 
Lemma  3.3.  Thus,  applying  Lemma  3.3  [with  ip  =  a+  and  X  =  B\],  we  get 


P  [(a+(Bi))*  >  x]  ~  ((r  -  7)  +  h'(oo))  >  *] 


(C.7) 


given  (4.1)  and  (C.6). 


(Step  2)  The  tail  of  asymptotics  of  M  will  be  identified  through  a  well-known  result  of  Veraverbeke 
[28,  Theorem  2(B),  p.  35].  To  prepare  for  it,  with  the  definition  (4.8)  of  A'j,  we  remark  the  inequality 
A’i  <  a(Bi)  +  6(0)  via  (C.3),  whence  A'+  <  a+(Bi)  +  b( 0)+,  and  the  integrability  of  A'+  is  implied 
by  that  of  a+(B\).  Also,  appealing  to  (C.3)  again,  we  get 

P  [A’i  >  .r]  =  P  [«+(Hi)-  |  b{h)  -  6(0)  |  +6(0)  >  *] 

~  P  [a+(Bi)  >  .r]  (C.8) 
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where  the  equivalence  is  validated  by  Lemma  3.1(2)  [with  A'  =  a+{B i),  Y  =|  6(0)  —  6(/f) 
d  =  6(0)]  under  the  independence  of  the  rvs  B\  and  I\ .  In  conclusion,  we  have 


and 


P  [a+  >  a-]  ~  P  [a+(Bi)  >  o’] 


(C.9) 


so  that  Xi  E  C  under  (H6).  Now,  the  same  arguments  based  on  Tchebychev’s  inequality  which 
gave  (C.6),  when  applied  to  the  rv  A'+,  yield 


0  <  E  [a+ 


<  OO. 


(C-10) 


Consequently,  the  integrated  tail  rv  (A+)*  associated  with  A+  is  well  defined,  and  the  equivalence 
(C.9)  gives 


roo  r  ..  ro o 

J  P  |A+  >  u  du  ~  J  P  [a+{Bi)  >  «]  du. 


(C.ll) 


Hence,  under  (C.6),  by  the  definition  of  the  integrated  tail  of  the  rv  a+(B\),  we  conclude  that 


ro o 

/  P  [A' i  >  u] 

J  X 


du 


roo 

/  P  [a+{Bi)  >  u]  du 

J  X 

E[a+(H!)]  P[(a+(B1))*>x] 


(C.12) 


~  ((r-  7)  +  h'(oo))  E  [B-t]  P  [a+(B^  >  *]  (C.13) 

where  the  asymptotic  equivalence  (C.13)  follows  from  (C.7). 

Using  (H7)  we  conclude  from  (C.13)  and  Lemma  3.1  that  the  rv  A+  has  an  integrated  tail  in  S. 
Since  E  [A'i]  <  0  under  (H5),  Theorem  2(B)  of  [28,  p.  35]  yields 

roo 

-E  [A'i]  P  [M  >x]~  P  [A'i  >  u]  du  (C.14) 

J  X 

and  upon  substituting  (C.13)  into  (C.14),  we  readily  obtain  the  asymptotics 

P  [M  >  x]  ~  K(h)P  [a+(B^)  >  x]  (C.15) 

where  K(h)  is  given  by  (4.11). 

(Step  3)  To  discover  the  tail  asymptotics  of  the  rv  Aq,  we  observe  from  (4.5)  that 

P  [A'0  >  x]  =  (1  -  p)P  [a(Bi)  +  b(Ii)  >  x]  +pP  [ a(B *)  +  6(0)  >|)  ,  x  >  0  (C.16) 

with  /*,  B\  and  B\  independent  rvs. 

Under  (H6),  the  rv  a+(B\)  belongs  to  C.  Thus,  recalling  (C.3)  and  applying  Lemma  3.1(2)  [with 
A  =  a+(Bi),  Y  =|  6(0)  —  6(/jf)  |  and  d  =  6(0)],  we  get 

P  [,/(/!,)  •  6(7; )  >  x  ~  P|,/-(/i|)  6(0)  6(  / 1* )  -6(0)  >.r; 

=  P  [«+(Bi)  >  .r]  .  (C.17) 
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Lemma  3.2  [with  A'  =  a+{B i)]  and  (C.7)  successively  yield 


P  [a+(5i)  >  x]  =  o(P  [(a+(Bi))*  >  x] ) 

=  o(P[a+(Bf)>  x]);  (C.18) 

the  required  moment  conditions  0  <  E  [Bi]  ,  E  [a+(_Bi)[  <  oo  in  Lemma  3.2  hold  owing  to  (4.1)  and 
(C.6).  It  is  now  plain  from  (C.17)  and  (C.18)  that 

P  [a(B i)  +  b(If)  >  x]  =  o( P  [a+  ( Bf )  >  .r]).  (C.19) 

On  the  other  hand,  a+{Bf )  belongs  to  5  under  (H7)  (thus  to  £),  whence  asymptotically  equivalent 
to  a+(Bf)  +  6(0),  i.e., 

P  [a(Bf)  +  6(0)  >  x]  ~  P  [a+(^)  >  x]  .  (C.20) 

Combining  (C.16),  (C.19)  and  (C.20)  we  find 

P  [A'0  >  x]  ~  pP  [a+(B^)  >  .r]  .  (C.21) 

(Step  4)  Collecting  (C.15)  and  (C.21),  we  readily  conclude  to  (4.10)  and  (4.11)  by  an  application 
of  Lemma  3.1(3)  [with  A'  =  M ,  T  =  A'o,  Z  =  a+(B\),  c\  =  K '(h),  c-j  =  1].  By  Lemma  3.1(1), 
membership  of  V  in  5  follows  from  that  of  a+{Bi)  in  S.  I 

D  A  proof  of  Proposition  6.3 

We  need  the  following  fact  later  in  the  proof: 

Lemma  D.l  Consider  an  1R -valued  rv  X  such  that  E  [A']  =  0  and 

E  [e0x]  <  oo,  |0|  <  6*  (D.l) 

for  some  6*  >  0.  Then,  there  exists  6 **  in  the  interval  (0 ,6*)  and  £  >  0  such  that 

E  [e0x]  <  :Se%  \8\  <  0 **.  (D.2) 

Proof:  Fix  6  in  1R..  It  is  a  simple  matter  to  check  the  identity 

edx  =  1  +  Ox  +  O'2  £  [  £  e0sd£  dt ,  x  G  R. 
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The  zero-mean  condition  implies 


E  [e*A']  =  1  +  62K{6)  with  K{6)  :=  E 

and  we  conclude 

E  [e0x]  <  1  +  62\K{6)\  <  e02\K[e)\. 

The  proof  of  (D.2)  is  now  completed  by  noting  that  condition  (D.l)  ensures  the  existence  of  some 
6 **  in  the  interval  (0,(9*)  such  that 


n 

1 

-K> 

^  CO 

CO 

CD 

Jo  ' 

Jo  J 

£  “ 


-  sup  \K{0)\  <  oo. 
2  |6»|<0** 


To  proceed  with  the  proof  of  Proposition  6.3,  we  set  C  :=  E  \C2\  >  0,  and  pick  scalars 

x  >  0,  a,  s  >  0,  T),  n  G  (0, 1)  and  0  G  (0,  C). 


We  also  write 

z(x) 

and  for  easy  reference,  we  set 

Mi  :  =  sup  Mi(6) 
o<e<e0 


(!  -v)„ 

- X 


and  n(x) 


I  (1  -  fl)z(x)  . 
1  C  +  a  J 


Mi(0 0)  and  W  :=  sup  Nt(6)  =  Nt(0o), 
0  <0<0O 


i  =  1,2. 


The  definition  of  the  rv  WA2'p2'  h  and  the  monotone  character  of  h  immediately  yield 

wA2,P2,-h  <  gup  (.4 2(Tn+1)  -  p2Tn  -  h(Tn ))  , 

4).  I  .... 

so  that  information  on  the  tail  of  the  rv  WA2,P2~h  can  in  principle  be  obtained  by  considering  the 
tail  of  the  maximum  associated  with  the  “perturbed  random  walk”  ^ .  1-2 (  1 )  —  p>  I t>  ~  h(  I „ )»  n  = 
0, 1, . . .}.  However,  the  matter  is  not  straightforward,  especially  when  L  =  0  in  (6.10),  as  we  note 
then  that 

lim  -  (A2(Tn+ 1)  -  p2Tn  -  h(Tn ))  =  0  a.s. 

n — ^00  72 

Although  a  Large  Deviations  Principle  will  hold  for  {A2(l„  .  1 )  —  p2Tn,  n  =  0, 1, . . .}  (and  even  for 
the  perturbed  random  walk  under  some  additional  conditions  on  h),  this  is  not  enough  to  guarantee 
exponential  decay  for  the  tail  of  the  rv  WA'2'p2~h.  The  basic  idea  of  the  proof  will  be  to  first 
“extract"  the  exponential  tails  associated  with  the  various  underlying  Large  Deviations  Principles. 
What  remains  will  provide  us  with  a  way  to  capture  the  effect  of  the  perturbation  h. 
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By  the  union  bound,  it  is  then  plain  that 

OO 

P  [ir'-'"-  >•  >  x]  <  £  P  [,12  (/•„..,)  -  p2Tn  -  h(T„)  >  x] 

n— 0 

<  A(x)  +  B(x)  (D.3) 

with 

OO 

A(x)  :=  Y  P  [A2{Tn+ 1)  -  (p2  +  e)Tn  >  i/x 

n— 0 

and 

OO 

B{x)  :=  ^2  P  [M{Tn+i}  -  p-iTn  -  h(Tn)  >  X ,  A2(Tn+i)  -  p2Tn  <  sTn  +  px]  . 

n— 0 

The  remainder  of  the  proof  consists  in  bounding  each  of  the  terms  A(x)  and  B(x). 


(Step  1)  Using  a  Chernoff  bound  argument  on  each  term  of  A(x)  and  the  independence  assumptions 
on  the  rvs  {(Cn,  T„),  n  =  1,2,...},  we  conclude  that 


A(x) 


OO 

p  [i i  >  Vx]  +  y  p 

n—  1 


}  n  ■  I  +  J^(h i  ~  (p2  +  ‘ )  O  )  >  Vx 

i=l 


< 


e-9r>x  Aq 


^1  +  N2 


J2  (E  [ee(y2-(P2  +  €)C2) 

n— 0 


o  <  e  <  e0. 


For  6  in  [0,#o],  r\{Q)  :=  E  [e0(V2  (p2+£)c'2)J  <  N2(6)  <  oo  under  (6.13),  and  we  get  ?’}(())  = 
— tE  [C2]  <  0  by  the  definition  of  p2.  Hence,  ri(6)  <  1  for  all  6  in  (0,#i)  with  0  <  6\  <  6q,  so  that 

A(x)  <  D^O)  e~6vx  with  D^O)  :=  i\q  ( 1  +  <  oo.  (D.4) 

V  1  -ri(0)J 


(Step  2)  To  handle  B(x)  we  note  the  following:  Whenever  A2(Tn+i)  —  p2Tn  —  h(Tn)  >  x  and 
A2{Tn+ 1)  -  p2Tn  <  -T„  px.  then  necessarily  A2{Tn+1)  -  p2Tn  >  h{Tn)  and  h(Tn)  +  x  <  sTn  +  px, 
the  latter  inequality  implying  z(x)  <  Tn+l-  Hence, 

OO 

B{x)  <  Y  P  [•  ^2 ( / ,, .  i )  -  p2l'„  >  h(Tn ),  z{x)  <  Tn+ 1] 

n— 0 

<  C(x)  +  D(x)  (D.5) 

with 

OO 

dx)  :=£p[(»  +  1  )(C  +  a)  +  pz{x)  <  Tn+ij 

n= 0 

and 

OO 

D(x)  :=  ^2  P  [A2{Tn+i)  -  p2Tn  >  //(  /'„).  z(x)  <  Tn+1  <  (n  +  1  ){C  +  a)  +  p,z(. r)]  . 

n=0 
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A  Chernoff  bound  argument  on  each  term  of  C{x)  yields 


C(x)  <  J2  e-eilz(x)V  |>(W+i-(n+i)(C'+«))j 

n= 0 

_  oo 

<  \I]C  Y,  r ■>(())" 

n— 0 

for  0  in  [0,  0o]  where  we  have  set  7-2(0)  :=  E  |^exp(@(C'2  —  C  —  a))j .  Since  ?’2 ( 0 )  =  0  and  r^O)  =  —a  < 

0,  there  exists  62  >  0  such  that  r 2(B)  <  1  on  (0,  62)-  Hence,  on  the  range  0  <  6  <  0^  :=  min{0i,  @2}, 
we  have 

H  f  —  0(C+a) 

C(x)  <  D2(0)  e~0lls^  with  D2(0)  :=  -p- - -w-  <  00.  (D.6) 

(1  -  r 2(0)) 

(Step  3)  Next,  we  note  that  the  condition  z(x)  <  Tn+ 1  <  (n  +  1)(C  +  a)  +  /. iz(x )  is  vacuous  unless 
(n  +  1)(C  +  a)  +  nz(x)  >  -.(x).  or  equivalently,  unless  n(x)  <  n,  whence  the  first  n(x)  terms  in  the 
sum  D(x)  equal  zero.  With  [3\  :=  C  —  /3  >  0,  we  get 

D(x)  <  E(x)  +  F(x)  (D.7) 

with 

OO 

E(x)  :=  £  P  lTn  <  n(h\ 

n=n(x) 

and 

OO 

F(x)  :=  Y  P  [A2(Tn+i)  -  p2Tn  >  h(Tn),n/3i  <  Tn\  . 

n=n(x) 

The  usual  Chernoff  bound  argument  now  yields 

OO 

E(x)  <  e6!h  r* W"-1  (D.8) 

n=n(x) 

for  all  0  >  0,  where  7-3(0)  :=  E  .  Since  7-3(0)  =  1  and  7-3(0)  =  —(3  <  0,  there  exists  @4  in 

(0,#o)  such  that  r%(0)  <  1  for  all  0  in  (0,  #4).  Hence,  on  the  range  0  <  0  <  0$  :=  min{03,04},  we 
deduce  from  (D.8)  that 

,  ,  ee'h 

E(x)<D3(0)r3(0)"t»l  with  D 3(0)  :=  (1  _  <  oo.  (D.9) 

(Step  4)  Collecting  the  bounds  (D.3),  (D.5)  and  (D.7),  we  find 

P  jir  *  ••/'•••  i'  >  3]  <  R(x)  +  F(x)  (D.10) 
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with 

R(x)  :=  A(x)  +  C(x)  +  E(x). 

Making  use  of  (D.4),  (D.6)  and  (D.9),  it  follows  that 

R(x)  <  1)00)  e~0r,x  +  D2(0)  e~6llz{x)  +  I):i(0)  r3(0)n ^ 

whenever  0  <  6  <  6 5.  Let  70  be  a  constant  such  that  n(x)  >  70X  for  x  large  enough,  e.g.,  select  70 
as  70  :=  (1  —  //)(1  —  rj)/{2s(C  +  a).  With  this  notation  we  get 

R(x)  <  D(0)  e~dx  (D.ll) 

for  x  large  enough,  where  0  <  d  :=  min  ( 61 7,  —70  \og(r3(0)),  0p(  1  —  i])/s)  and  D{6)  :=  Yl'i-i  F>i(0). 

(Step  5)  For  the  last  term  F(x ),  observe  from  the  monotone  character  of  h  that  h(Tn)  >  h(nf3\)  if 
n/3 1  <  whence 

OO 

F(x)<  £  P  [A2(Tn+1)  -  p2Tn  -  h(ndi)  >  0\  .  (D.12) 

n=n(x) 

Pick  0  in  the  interval  (0,  6q).  For  each  n  =  1,  2, . . .,  a  Chernoff  bound  argument  gives 
P  [A2(Tn+i)  -  p2Tn  -  h(n/3i)  >  0] 

<  g-dh(np i)jj  j’gfi'i'n+ij  (>■>.<  9.)'  "  1 

<  A 1  .V2c  (D.13) 

By  Lemma  D.l,  under  the  moment  relation  E  [I2]  =  p2E[Cf2],  there  exist  £  >  0  and  0q  in  (0,  #o) 
such  that 

K  [r0(V"  <  (A°\  0  (r  (0. //,;).  (D.14) 

Combining  (D.13)  and  (D.14),  we  find 

P  [A2(Tn+ 1)  -  p2Tn  -  //(//hi  >  0]  <  JVj^e-^^T02).  (D.15) 

This  last  upper  bound  is  best,  i.e. ,  smallest,  for  0  =  h(/3in)/n£,  a  quantity  that  can  be  made 
arbitrary  small  since  (6.10)  holds  with  L  =  0,  hence  smaller  than  0q,  for  n  large  enough,  in  which 
case  (D.15)  can  be  tightened  to 

h(nf)  |  )2 

P  [A2(Tn+ 1)  -  P2Tn  -  1101/30  >  0]  <  N3N2e  2<  .  (D.16) 

In  particular,  for  x  large  enough,  we  conclude  from  (D.12)  and  (D.16)  that 

“  feO/3,)2 

/•'(./•)  <  A"  1  Ah  e  (D.17) 

n=n(x) 
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Consequently,  for  x  large  enough,  we  have  n(x)  >  7o:r,  and  with  71  :=  /4i7o  and  72  :=  /3 i/2£,  a 
standard  bounding  argument  shows  that 


y  e 

n=n(x) 


<  E  * 

n=n(x) 


729(nPl) 
n/3 1 


< 


L 


-I2g(pit)2 / (ht  ^-1- 


70  a? 


kl 


jix 


e-72  9(t?lt(U  .  _L7>>  , 

Pi 


(D.18) 


as  we  note  that  f  — *■  #(f)/ v7  =  infs>t  /i(s)/ 7*  is  nondecreasing.  Reporting  (D.18)  into  (D.17)  yields 
F(x)  <  MRh(x)  with  M  :=  N1N2/P1,  so  that  (D.10)  and  (D.ll)  together  now  imply 

P  h  >  x\  <  D(6)  e~dx  +  MRg(x )  (D.19) 

for  x  large  enough. 


(Step  6)  We  are  now  in  position  to  prove  (6.14)  and  (6.15).  In  view  of  (D.19),  (6.14)  will  fol¬ 
low  with  73  =  M  if  we  show  that  lim.,.  .  N  t  ttX/Rg(.r)  =  0  for  any  a  >  0,  or  equivalently  that 
limaMOO  Rg(x)eax  =  oo  for  any  a  >  0.  To  that  end,  with 

\  0,  q(t )2 

G(t)  := — t  72 7  1  *  >  0 

7i  t 

we  get 

poo 

Rg(x)eax  =  eax  eG^e~at^dt 
J  ■y  l  x 

>  (  inf  •  f°°  eax~atl‘ 'n  dt  =  —  inf  eG'(t). 

\t> -fix  )  J~nx  a  *>71® 


The  desired  conclusion  is  obtained  if  we  show  that  Imp  ,rx  G{t)  =  too,  a  fact  which  follows  from 
L  =  0  in  (6.10),  the  fact  g(t)/t  <  h(t)/t  ( t  >  0)  and  the  identity 


G(t)=t(~-  72 
\7i 


t  >  0. 


We  now  turn  to  the  proof  of  (6.15).  For  any  given  constant  C  >  0,  assumption  (2)  implies  that 
h(t)2 /t  >  C't2u°  >  Ctv  for  t  large  enough  whenever  v  lies  in  (0,  2vq).  Therefore,  for  x  large  enough, 
a  simple  bounding  argument  and  the  well-known  asymptotics  of  the  incomplete  Gamma  function 
lead  to 

poo 

Rg(x)  <  /  e~sr  dt  ~  60  (71ar)1"*'  (D.20) 

J-yix 

with  6  :=  C72  and  :=  Combining  (6.14)  and  (D.20)  gives  (6.15)  with  6\  =  67^,62  = 

73<5o7i1_'y)-  I 
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E  A  proof  of  Proposition  8.1 


We  show  that  the  choices  h  =  U  and  h  =  —  U  meet  the  requirements  of  Claims  1  and  2,  respectively. 
To  avoid  unnecessary  repetitions,  we  start  with  some  comments  that  are  common  to  both  Claims. 

(Fact  1):  Because  lim,.  .-x:  ?»/;  (.r)  =  oo  by  Lemma  3.2,  the  same  property  holds  for  u  by  the 
asymptotic  equivalence  ms1  ~  u,  hence  for  U.  The  mapping  U  being  concave  and  increasing  on 
Dl  .  .  the  convergence 

lim  LSll  =  ijm  =  lim  u1  (t)  =:  v!(oo)  (E.l) 

t — >oo  t  t — >oo  t  t — >oo 

takes  place  in  a  monotonically  decreasing  manner,  whence  the  limit  exists  and  is  finite. 

(Fact  2):  With  (5.2)  and  (5.4)  in  mind  for  the  situations  at  hand,  for  every  ij  in  1R.  we  write 

av(x)  :=  ax  +  7jU(x}r  x  >  0,  (E.2) 

and 

A'?  :=  aBx  -  (c  -  p2)h  +  +  U(h )).  (E.3) 

By  Jensen’s  inequality  we  have  0  <  E[C(i?i)]  <  U(E[Bi])  <  oo,  with  similar  inequalities  for  Ji, 
whence  the  rvs  U{B\)  and  U{I\)  are  both  integrable.  Thus,  the  rvs  av(Bi)  and  X'l  are  integrable; 
we  have  E  [A'j?]  <  0  for  small  enough  //  in  view  of  the  fact  that  lirn^jo  E  [A')7]  <  0  by  the  stability 
condition. 

(Fact  3):  For  a  given  p  in  II,  whenever  the  mapping  av  is  strictly  increasing  in  the  limit ,  we  have 

P  [a+(Bi)  >  x\ 
lim  — — — - rT 

x— oo  p  [aBi  >  x\ 


by  virtue  of  Lemma  3.6  since  B\  G  MDA(  A)  (and  ms1  ~  U).  A  similar  argument,  making  use  this 
time  of  the  consequence  (3.12)  of  Lemma  3.6,  also  yields 


P  a+(Bi)>a+(x) 


lim  r 

^°°  P  [cvBi  >  a+(x) 


P  [Bi  >  x] 


lim  r 

P  [aB1  >  an(x) 


=  e 


(E.4) 


P  a^(B^)  >  x 

lim  — 

P  [aBi  >  x\ 


lim  — 

x—>oo  p 


P  [B\  >  $] 


aB\  >  a] \  (x) 


=  e 


(E.5) 


In  other  words,  whenever  the  mapping  at]  is  strictly  increasing  in  the  limit,  we  conclude  from  (E.4) 
that  the  rvs  aB\  and  a+(B\)  are  both  tail  equivalent,  while  (E.5)  implies  that  the  rvs  aB\  and 
a+(Bi)  are  tail  equivalent.  By  Lemma  3.1  (1)  the  rvs  a+{B\)  and  a+{B\)  are  elements  of  C  and 
5,  respectively. 

We  are  now  ready  to  discuss  Claims  1  and  2. 
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(Claim  1):  Define  h  :  =  U  and  fix  e  >  0.  The  mapping  h  is  subadditive  and  absolutely  continuous 
by  concavity  (with  h( 0)  =  0).  The  scaled  mapping  eh  automatically  satisfies  (H3)  and  (H4)  (under 
the  second  half  of  (8.1)).  From  this  latter  condition  we  conclude  via  Fact  3  that  af  {B i)  and  af  (Bf) 
are  elements  of  C  and  5,  respectively,  i.e.,  conditions  (H6)  and  (H7)  both  hold  for  eh.  Conditions 
(HI),  (H2)  and  (H3)  on  eh  are  equivalent  to 

c  —  (?’i  +  p-j)  <  eU'(x)  <  c  —  P2  a.e.  on  1R.+  (E.6) 

with  U'(oo)  finite.  Because  U  is  a  non-decreasing  and  concave  function,  it  follows  that  0  <  U'( 00)  < 
U'(x)  <  U'( 0+)  a.e.  on  IR+  with  U'( 0+)  finite  by  Lemma  8.1.  The  constraints  (E.6)  are  therefore 
implied  by  requiring 

c  —  (?’i  +  P2)  <  eU'(oo)  <  eU'(0+)  <  c  —  p2  a.e.  on  R+  (E.7) 

and  under  (8.1)  this  is  obviously  satisfied  if  r  is  chosen  sufficiently  small.  If  t  is  taken  small  enough, 
we  see  from  the  discussion  above  that  (H5)  holds  as  well. 

(Claim  2):  This  time,  define  h  :=  —U  and  fix  e  >  0.  The  absolutely  continuous  mapping  h  is  now 
superadditive.  Since  U  is  non-decreasing,  the  scaled  mapping  eh  automatically  satisfies  (H2),  while 
condition  (H3)  holds  by  the  remarks  leading  to  (E.l).  Moreover,  by  Fact  2  we  see  that  eh  satisfies 
(Ho). 

The  mapping  a_e  will  be  strictly  increasing  if  a'_£{x)  >  0  a.e.  on  R+,  or  equivalently, 

a'_s{x )  =  r  1  +  P2  —  c  —  sU'(x)  >  0  a.e.  on  ]R.+  .  (E.8) 

By  remarks  made  in  the  proof  of  Claim  1,  this  last  requirement  will  hold  if 

a1  ,.(x)  >  ?’i  +  p2  —  c  —  e'Lr,(0+)  >  0  a.e.  on  ]R+  (E.9) 

by  virtue  of  the  concavity  of  U.  This  is  always  possible  owing  to  (8.1)  by  selecting  e  sufficiently 
small,  in  which  case  a-s  is  strictly  increasing  on  ]R+  with  lim  ,.  ..x.  a-£{x)  =  00. 

Consec|uently,  for  i  sufficiently  small,  it  is  the  case  that  eh  satisfies  (H4),  hence  (H6)  and  (H7)  by 
appealing  to  Fact  3.  I 
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